Astronomy & Astrophysics manuscript no. m09bnl0 © ESO 2013 

April 16, 2013 



Chemistry of massive young stellar objects with a disk-like 

structure 

K. Isokoski 1 , S. Bottinelli 23 , and E. F. van Dishoeck 45 

'Raymond and Beverly Sackler Laboratory for Astrophysics, Leiden Observatory, Leiden University, PO Box 95 13, 2300 RA Leiden, 
The Netherlands 

2 Universite de Toulouse, UPS-OMP, IRAP, Toulouse, France 

3 CNRS, IRAP, 9 Av. colonel Roche, BP 44346, F-31028 Toulouse cedex 4, France 

4 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 

5 Max-Planck Institut fur Extraterrestrische Physik (MPE), Giessenbachstr. 1, 85748 Garching, Germany 



m 



ABSTRACT 



in 



43 
Oh 



f— ^ Preprint online version: April 16, 2013 

(N 

Sh 

Aims. Our goal is to take an inventory of complex molecules in three well-known high-mass protostars for which disks or toroids 

■^s| have been claimed and to study the similarities and differences with a sample of massive YSOs without evidence of such flattened 

disk-like structures. With a disk-like geometry, UV radiation can escape more readily and potentially affect the ice and gas chemistry 
on hot-core scales. 
Methods. A partial submillimeter line survey, targeting CH 3 OH, H 2 CO, C 2 H 5 OH, HCOOCH 3 , CH 3 OCH 3 , CH 3 CN, HNCO, 

, , NH 2 CHO, C 2 H 5 CN, CH 2 CO, HCOOH, CH 3 CHO, and CH 3 CCH, was made toward three massive YSOs with disk-like structures, 

^ IRAS20 126+4 104, IRAS18089-1732, and G31.41+0.31. Rotation temperatures and column densities were determined by the rotation 

diagram method, as well as by independent spectral modeling. The molecular abundances were compared with previous observations 

\^ of massive YSOs without evidence of any disk structure, targeting the same molecules with the same settings and using the same 

analysis method. 

Results. Consistent with previous studies, different complex organic species have different characteristic rotation temperatures and 
can be classified either as warm (>100 K) or cold (<100 K). The excitation temperatures and abundance ratios are similar from source 
to source and no significant difference can be established between the two source types. Acetone, CH 3 COCH 3 , is detected for the first 
time in G31.41+0.31 and IRAS 18089- 1732. Temperatures and abundances derived from the two analysis methods generally agree 

"£i within factors of a few. 

Conclusions. The lack of chemical differentiation between massive YSOs with and without observed disks suggest either that the 

1 1 chemical complexity is already fully established in the ices in the cold prestellar phase or that the material experiences similar physi- 

cal conditions and UV exposure through outflow cavities during the short embedded lifetime. 
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£^ 1. Introduction different beams, different frequency ranges and different analysis 

' techniques. 
Tj" Millimeter lines from complex organic molecules are widely as- 

O sociated with high-mass star forming regions and indeed form Chemical abundances depend on the physical structure of 
Cn one of the signposts of the deeply embedded phase of star for- the source such as temperature, density and their evolution with 
7"J mation (e.g., Blake et al. 1987; Hatchell et al. 1998; Gibb et al. time, as well as the amount of UV radiation impinging on the gas 
j> 2000; Fontani et al. 2007; Requena-Torres et al. 2008; Belloche and dust. In contrast with the case of solar-mass stars, the phys- 
et al. 2009; Zernickel et al. 2012). Many studies of the chem- ical structures and mechanisms for forming massive (M > 8 
/\ istry in such regions have been carried out, either through com- M ) stars are still poorly understood. Indeed, theoretically, the 
H plete spectral surveys of individual sources or by targetting in- powerful UV radiation pressure from a high-mass protostellar 
dividual molecules in a larger number of sources (see Herbst & object (HMPO) should prevent further accretion and so inhibit 
van Dishoeck 2009; Caselli & Ceccarelli 2012, for reviews). In the formation of more massive stars (Zinnecker & Yorke 2007). 
spite of all this work, only few systematic studies of the abun- However, a number of recent studies have claimed the presence 
dances of commonly observed complex molecules have been of disk- or toroid-like 'equatorial' structures surrounding a hand- 
performed across a sample of massive YSOs, to search for sim- ful of high-mass protostars (Cesaroni et al. 2007). These data 
ilarities or differences depending on physical structure and evo- support theories in which high-mass star formation proceeds in a 
lutionary state of the object. Intercomparison of published data similar way as that of low-mass stars via a disk accretion phase in 
sets is often complicated by the use of different telescopes with which high accretion rates and non spherically symmetric struc- 
tures overcome the problem of radiation pressure. The best ev- 

idence so far is that for a ~5000 AU disk in Keplerian rotation 

Send offprint requests to: K. Isokoski; e-mail: around IRAS20126+4104, claimed by Cesaroni et al. (1999) on 

isokoski@strw.leidenuniv.nl the basis of the presence of a velocity gradient in the CH3CN 



Isokoski et al.: m09bnl0 




Disk 




Fig. 1: Illustration of a protostar with a spherical structure (left) 
and a protostar with a flattened disk-like structure (right) with 
enhanced UV photons illuminating the walls of the outflow cav- 
ity. 



emission perpendicular to the direction of the outflow, as pre- 
dicted by the disk-accretion paradigm. Surprisingly, even in the 
best case of IRAS20126+4104, the detailed chemistry of these 
(potential) disks has not yet been studied. 

Most chemical models invoke grain surface chemistry to cre- 
ate different generations of complex organic molecules (Tielens 
& Charnley 1997). Hydrogenation of solid O, C, N and CO dur- 
ing the cold (Td < 20 K) prestellar phase leads to ample pro- 
duction of CH3OH and other hydrogenated species like H2O, 
CH4 and NH 3 (Tielens & Hagen 1982). Exposure to UV radi- 
ation results in photodissociation of these simple ices, with the 
fragments becoming mobile as the cloud core heats up during 
the protostellar phase. First generation complex molecules re- 
sult from the subsequent recombination of the photofragments, 
and will eventually evaporate once the grain temperature rises 
above the ice sublimation temperature of ~100 K (Garrod & 
Herbst 2006; Garrod et al. 2008). Good examples are C 2 H 5 OH, 
HCOOCH3 and CH3OCH3 resulting from mild UV processing 
of CH3OH ice (Oberg et al. 2009). Finally, a hot core gas phase 
chemistry between evaporated molecules can drive further com- 
plexity in second generation species (e.g., Millar et al. 1991; 
Charnley et al. 1992, 1995). 

One of the most obvious consequences of an equatorial 
rather than spherical structure is that UV radiation can more eas- 
ily escape the central source and illuminate the surface layers of 
the surrounding disk or toroid as well as the larger scale envelope 
(Bruderer et al. 2009, 2010) (Fig. 1). This can trigger enhanced 
formation of complex organic molecules in the ices relative to 
methanol. Another effect of UV radiation is that it drives in- 
creased photodissociation of gaseous N2 and CO. The resulting 
atomic N and C would then be available for grain surface chem- 
istry potentially leading to enhanced abundances of species like 
HNCO and NH 2 CHO. 

To investigate the effects of disk-like structures on the chem- 
istry, we present here a single-dish survey using the James Clerk 
Maxwell Telescope of three HMPOs for which large equato- 
rial structures (size> 2000 AU) have been inferred, namely 
IRAS20126+4104, IRAS 18089-1732, and G31.41+0.31. The 
results are compared with those of a recent survey of a sam- 
ple of HMPOs by Bisschop et al. (2007) (hereafter BIS07), tar- 
geting many of the same molecules and settings. The use of 
the same telescope and analysis method allows a meaningful 
comparison between the two samples of sources. BIS07 found 
that the O-rich complex molecules are closely correlated with 
the grain surface product CH3OH supporting the above general 
chemical scenario. N-rich organic molecules do not appear to 
be correlated with O-rich ones, but overall, the relative abun- 



dances of the various species are found to be remarkably con- 
stant within one chemical family. One of the main questions is 
whether this similarity in abundance ratios also holds for sources 
with disk-like structures. Although our observations do not spa- 
tially resolve these structures, they are sensitive enough and span 
a large enough energy range to determine abundances on scales 
of ~l" and thereby set the scene for future high-angular reso- 
lution observations with interferometers like the Atacama Large 
Millimeter/submillimeter Array (ALMA). Moreover, current in- 
terferometers resolve out part of the emission, which is why 
single-dish observations are still meaningful. 

This paper is organized as follows. In §2, the observed 
sources and frequency settings are introduced and the details 
of the observations are presented. §3 focuses on the data anal- 
ysis methods. Specifically, two techniques are used to determine 
excitation temperatures and column densities: the rotation dia- 
gram method employed by BIS07 and spectral modeling tools 
in which the observed spectra are simulated directly. §4 presents 
the results from both analysis methods and discusses their advan- 
tages and disadvantages. Finally, §5 compares our results with 
those of BIS07 to draw conclusions on similarities and differ- 
ences in chemical abundances between sources with and without 
large equatorial structures. 

2. Observations 

2.1. Observed sources 

Table 1 gives the coordinates, luminosity L, distance d, galacto- 
centric radius Rgc, velocity of the local standard of rest Vlsr 
and the typical line width for the observed sources. The se- 
lected sources are massive young stellar objects (YSOs), for 
which strong evidence exists for a circumstellar disk structure. 
All sources are expected to harbor a hot core: a compact, dense 
(>10 7 cm 4 ) and warm (>100 K) region with complex chem- 
istry triggered by the grain mantle evaporation (Kurtz et al. Kurtz 
et al. (2000)). CH3CN emission from >100 K gas is present in 
all sources. Moreover, CH3OH 7^-6^ transitions (338.5 GHz) 
with main beam temperatures of > 1 K have been observed for 
these sources. Sources also needed to be visible from the James 
Clerk Maxwell Telescope (JCMT) ' . 

2.1.1. IRAS201 26+41 04 

IRAS20126+4104 (hereafter IRAS20126) is a luminous (~10 4 
L Q ) YSO located relatively nearby at a distance of 1.64 kpc 
(Moscadelli et al. 201 1). It was first identified in the IRAS point 
source catalog by the IR colours typical of ultracompact HII re- 
gions and by H2O maser emission characteristic of high-mass 
star formation (Comoretto et al. 1990). IRAS20126 features a 
~0.25-pc scale inner jet traced by H2O maser spots in the SE- 
NW direction with decreasing velocity gradient (Tofani et al. 
1995). Source and masers are embedded inside a dense, parsec- 
scale molecular clump (Estalella et al. 1993; Cesaroni et al. 
1999). The inner jet feeds into a larger scale bipolar outflow 
with the two having reversed velocity lobes (Wilking et al. 1990; 
Cesaroni et al. 1999). The reversal is likely to be due to preces- 
sion of the inner jet caused by a companion separated by a dis- 
tance of -0.5" (850 AU) (Hofner et al. 1999; Shepherd et al. 



The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre, on behalf of the Particle Physics and Astronomy 
Research council of the United Kindom, the Netherlands Organization 
for Scientific Research and the National Research Council of Canada. 
The project ID is m09bnl0. 
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Table 1 : Coordinates, luminosity, distance, galactocentric radius, velocity of the local standard of rest, 12 C/ 13 C isotope ratio and the 
typical line width for the observed sources. 



Sources 


a(2000) 


(5(2000) 


L 


d 


^GC* 


Vlsr 


12 C/"C 


AV 








[10 5 L e ] 


[kpc] 


[kpc] 


[kms- 1 ] 




[kms- 1 ] 


IRAS20126+4104 


20:14:26.04 


+41:13:32.5 


0.13" 


1.64" 


8.3 


-3.5 


70 


6 C 


IRAS18089-1732 


18:11:51.40 


-17:31:28.5 


0.32 rf 


2.34 e 


6.2 


33.8 


54 


5 C 


G31.41+0.31 


18:47:34.33 


-01:12:46.5 


2.6' 


7.9 s 


4.5 


97.0 


41 


6-10'" 



(a) Cesaroni et al. (1997), (b) Moscadelli et al. (201 1), (c) Leurini et al. (2007), (d) Sridharan et al. (2002), (e) Xu et al. (201 1), (f) Cesaroni et al. 
(1994b), (g) Churchwell et al. (1990), (h) Fontani et al. (2007). 

* The galactocentric radii were calculated using distances d in this table and a IAU recommended distance from the galactic center R = 8.5 kpc. 
The 12 C/' 3 C isotope ratios are calculated from Eq. 1 1 . 



2000; Cesaroni et al. 2005; Sridharan et al. 2005). A rotating, 
flattened, Keplerian disk structure has been detected perpendic- 
ular to the inner jet. Observations of CH3CN transitions show a 
Keplerian circumstellar disk (radius ~1000 AU) with a velocity 
gradient perpendicular to the jet and a hot core with a diameter 
of ~0.0082 pc and a temperature of ~200 K at a geometric center 
of the outflow (Cesaroni et al. 1997; Zhang et al. 1998; Cesaroni 
et al. 1999). Direct near-infrared (NIR) observations show a disk 
structure as a dark line (Sridharan et al. 2005). The disk shows 
a temperature and density gradient and is going through infall 
of material with a rate of ~2 x 10~ 3 M Q yr -1 as expected for 
a protostar of this mass and luminosity (Cesaroni et al. 2005). 
Recent modeling of the Spectral Energy Distribution (SED) of 
IRAS20126 shows indeed a better fit when a disk is included 
(Johnston et al. 2011). 



2.1.2. IRAS18089-1732 

IRAS 18089-1732 (hereafter IRAS18089) is a luminous 
(~10 45 L Q , Sridharan et al. 2002) YSO located at a distance 
of 2.34 kpc (Xu et al. 2011). It was identified based on CS 
detections of bright IRAS point sources with colours similar 
to ultracompact HII regions and the absence of significant 
free-free emission (Sridharan et al. 2002), and with H2O and 
CH3OH maser emission (Beuther et al. 2002). The CO line 
profile shows a wing structure characteristic of an outflow, 
although no clear outflow structure could be resolved from the 
CO maps (Beuther et al. 2002). A collimated outflow in the 
Northern direction is however seen in SiO emission on scales of 
5" (Beuther et al. 2004). IRAS 18089 has a highly non-circular 
dust core of ~3000 AU diameter (~1"), with optically thick 
CH 3 CN at -350 K (Beuther et al. 2005). HCOOCH3 was 
found to be optically thin, with emission confined to the core, 
and showing a velocity gradient perpendicular to the outflow 
indicative of a rotating disk (Beuther et al. 2004). Also hot 
NH3 shows a velocity gradient perpendicular to the outflow, 
although no Keplerian rotation was found, possibly due to 
infall and/or self gravitation (Beuther & Walsh 2008). Several 
hot-core molecules (HCOOCH3, CH3CN, CH3OCH3, HNCO, 
NH 2 CHO, CH3OH, C 2 H 5 OH) were mapped by Beuther et al. 
(2005) but no column densities or abundances were reported. 



2.1.3. G31.41+0.31 

G31.41+0.31 (hereafter G31) is a luminous (2.6 x 10 5 L G , 
Cesaroni et al. 1994b) YSO at a distance of 7.9 kpc (Churchwell 
et al. 1990). Preliminary evidence for a rotating disk with a 
perpendicular bipolar outflow was reported by Cesaroni et al. 
(1994a,b) and Olmi et al. (1996) showing a velocity gradient 



across the core in the NE-SW direction, similar to previously 
detected OH masers (Gaume & Mutel 1987). High-angular res- 
olution CH3CN observations by Beltran et al. (2005) could not 
detect Keplerian rotation typical for less luminous stars, never- 
theless a toroidal structure undergoing gravitational collapse and 
fast accretion (~3 x 10~ 2 M Q yr -1 ) onto the central object was 
found. The G31 hot molecular core (HMC) is part of a complex 
region in which multiple stellar sources are detected (Benjamin 
et al. 2003); indeed, it is separated from an ultracompact (UC) 
HII region by only ~5" and overlaps with a diffuse halo of free- 
free emission, possibly associated with the UC HII region itself 
(Cesaroni et al. 1998). More recent interferometric observations 
confirm the velocity gradient in the NH3 (4,4) inversion transi- 
tion and in CH3CN data (Cesaroni et al. 2010, 201 1). Line pro- 
files look like a rotating toroid with infall motion. Several com- 
plex hot-core molecules have been observed in G31, including 
glycolaldehyde CH 2 OHCHO (Beltran et al. 2009), but again no 
abundance ratios have been presented. 

2.2. Observational details 

The observations were performed at the JCMT on Mauna Kea, 
Hawaii, between August 2007 and September 2009. The obser- 
vations of the 338 GHz region covering CH3OH (Ik — > 6a:) 
transitions were taken from JCMT archive. The front ends con- 
sisted of the facility receivers A3 (230 GHz region) and HARP- 
B (340 GHz region). The back-end was the digital autocorrela- 
tion spectrometer (ACSIS), covering 400 and 250 MHz of in- 
stantaneous bandwidth for A3 and HARP-B, respectively, with 
a channel width of 50 kHz. The noise level for both receivers 
was r rms ~20 mK on a T* K scale when binned to 0.5 km s . The 
integration time was ~1 hr and 1.8 hr for A3 and HARP-B, re- 
spectively. The spectra were scaled from the observed antenna 
temperature, 7^, to main-beam temperature, Tmb, using main 
beam efficiencies 77MB of 0.69 and 0.63 at 230 GHz and 345 GHz, 
respectively. We adopt a T* A calibration error of 20%. 

The HPBW (half -power beam width, 6> B ) for the 230 and 
345 GHz band observations are 20-21" and 14", respectively. 
Emission from a volume with a source diameter 6s <6% under- 
goes beam dilution described by the beam-filling factor, tjsf'. 



«l 



°S +6 B 



(1) 



Table 2 gives the observed frequency settings and the tar- 
geted molecular lines. The settings were taken from BIS07 and 
were chosen to cover at least one strong line for the target 
molecules as well as lines of other interesting species. Strong 
lines of target molecules were chosen due to their high main- 
beam temperatures and minimum line-confusion in line surveys 
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of Orion-KL by Sutton et al. (1985) and Schilke et al. (1997) at 
230 GHz and 345 GHz, respectively. In order to determine rota- 
tion temperatures, we covered at least two transitions for a given 
species with E up <\00 K and £ up >100 K each. BIS07 used the 
single pixel receiver B3, the predecessor of HARP-B, together 
with DAS (Digital Autocorrelation Spectrometer) as the back- 
end, covering a larger instantaneous bandwidth of 500 MHz. 
Only the central receptor of HARP-B array is analysed here as 
no significant off source emission was detected in the complex 
molecules. 

Data reduction and line fitting were done using the CLASS 
software package 2 . Line assignments were done by compari- 
son of observed frequencies corrected for source velocity with 
the JPL 3 , CDMS 4 and NIST 5 catalogs (Pickett et al. 1998; 
Miiller et al. 2005). The line assignment/detection was based on 
Gaussian fitting with the following criteria: (i) the fitted line po- 
sition had to be within ±1 MHz of the catalog frequency, (ii) 
the FWHM had to be consistent with those given in Table 1 and 
(hi) a S IN >3 is required on the peak intensity. Sect. A.l in 
the Appendix describes in more detail the error analysis. In gen- 
eral, our errors on the integrated intensities are conservative and 
suggest a lower S /N than that on the peak intensity or obtained 
using more traditional error estimates. 



2 CLASS is part of the GILDAS software package developed by 
IRAM. 

3 http : //spec . jpl . nasa . gov/f tp/pub/catalog/catf orm . 
html 

4 http://www.phl .uni-koeln.de/vorhersagen 

5 http : //physics .nist . gov/PhysRefData/Micro/Html/ 
contents.html 



Isokoski et al.: m09bnl0 
Table 2: Observed frequency settings and molecular lines. 



Molecule 


Freq. 




fs 


Transition 


Additional molecules 




[GHz] 


[K] 


[D 2 ] 






CH 3 CN 


239.1195 


144.77 


811.86 


13* -12* 


CH 3 1J CN, HCOOCH3, CH3OCH3 


HNCO 


219.7983 


58.02 


28.112 


10o.10.Il -9o,9,10 


H 2 13 CO,C 2 H 5 CN 




240.8809 


112.53 


30.431 


Hi, 11, 12 - 10i, 10,11 


CH3OCH3, CH3OH, HN 13 CO 


HCOOCH3 


222.3453 


37.89 


42.100 


85,4 - 7 4 3 


CH3OCH3, NH 2 CHO 




225.2568 


125.50 


33.070 


186,12 - 176.12 


H 2 CO, CH3OCH3, 13 CH30H 


H 2 CO 


364.2752 


158.42 


52.165 


5 3 .3 - 4 3 ,2 


C2H5OH 


CH3CN 


331.0716 


151.11 


513.924 


18* -17* 


HCOOCH,, HNCO, CHi 13 CN 


HNCO 


352.8979 


187.25 


43.387 


16l, 15,17 - 15l,14,16 


C2H5CN, C2H5OH, HCOOCH3 


NHoCHO 


345.1826 


151.59 


664.219 


17o,17 - 16o,16 


HCOOCH3, C,H,OH, 13 CH30H 


HCOOCH3 


354.6084 


293.39 


87.321 


332,32 - 322,31 


NH 2 CHO, C2H5CN 



3. Data analysis 

3.1. Rotation diagrams 

Rotation temperatures and column densities were obtained 
through the rotation diagram (RTD) method (Goldsmith & 
Langer 1999), when 3 or more lines are detected over a suffi- 
ciently large energy range. Integrated main-beam temperatures, 
f TmbcIV, can be related to the column density in the upper en- 
ergy level by: 



JV up _ 3kjT MB dV 
g up 8n 3 vfi 2 S 



(2) 



where N up is the column density in the upper level, g up is the 
degeneracy of the upper level, k is Boltzmann's constant, v is the 
transition frequency, fi is the dipole moment and S is the line 
strength. The total beam-averaged column density Nt in cm -2 
can then be computed from: 



Nup 

gup 



N 7 



Q(J rot ) 



£• up / ' 10 



(3) 



where Q(T rot ) is the rotational partition function, and E np is the 
upper level energy in K. 

Blended transitions of a given species with similar 
E up (A£ , up <30 K) were assigned intensities according to their 
Einstein coefficients (A) and upper level degeneracies {g up ): 



jT MS dV(i)= f: 



T MB dV x 



A'gi 



up 



Z/A'sip 



where the summation is over all the contributing transitions. 
Blended transitions with different E up , or with contamination 
from transitions belonging to other species, were excluded from 
the RTD fit. 

The beam averaged column density, Nj, is converted to the 
source-averaged column density Ns using the beam-filling fac- 
tor 77 BF : 



N s 



N T 



(5) 



The standard RTD method assumes that the lines are opti- 
cally thin. Lines with strong optical depth, determined from the 
arguments in §4.1 as well as the models discussed in §3.2, were 
excluded from the fit. For CH3OH, also low-£up lines arising 
from a cold extended component (see Sect. 4.4. 1) were excluded. 



Differential beam dilution is taken into account by multi- 
plying the line intensities in the 230 GHz range by ^ BF (340 
GHz)/77bf(230 GHz). For warm and cold molecules, beam di- 
lution is derived assuming source diameters 6't=iook (see Eq. 10 
below) and 14", respectively (see Eq. 1). All emission is thus 
assumed to be contained within the smallest beam size. Same 
approach was used in BIS07. 

The vibrational partition function was ignored assuming that 
all emission arises from the ground vibrational state so that the 
vibrational partition function can be set to unity. This approxi- 
mation can cause an underestimation of the derived column den- 
sities for larger molecules, even at temperatures of 100-200 K. 
Indeed, Widicus Weaver et al. (2005) show that the error can 
be up to a factor of 2 for temperatures up to 300 K. Since this 
approximation affects all complex molecules (albeit at different 
levels), the overall error in abundance ratios will be less than a 
factor of 2 and well within the other uncertainties. 

3.2. Spectral modeling 

The alternative method for analysing the emission is to model 
the observed spectra directly. For this purpose, we used the so- 
called 'Weeds' extension of the CLASS software package , de- 
veloped to facilitate the analysis of large millimeter and sub- 
millimeter spectral surveys (Maret et al. 2011). In this model, 
the excitation of a species is assumed to be in LTE (Local 
Thermodynamic Equilibrium) at a temperature r ex . The bright- 
ness temperature, 7b, of a given species as a function of the rest 
frequency v is then given by: 



(4) T B (v) = 77BF [Jv (T ra ) - Jy (r bg )] ( 1 - e-™) , 



(6) 



where t]bf is the beam filling factor (see Eq. 1), J v is the radiation 
field such that: 



Jy(T) 



hv/k 



e hv/kT _ l 



and 7\,g is the temperature of the background emission. The 
HPBW 6b is calculated within the Weeds model as 6b = 
1.22c/vD, where c is the speed of light and D is the diameter 
of the telescope 7 . The opacity t(v) is: 



T(V) 



Nt 



8^v 2 g(r ex ) 



I>^P' 



e -K e /kT c 



w/^-l)^' 



(7) 



CLASS is part of the GILDAS software package developed by 
IRAM. 

7 For JCMT with 15-m antenna diameter, the equation gives a 6* B of 
21.9" and 14.8" for the 230 GHz band and 345 GHz bands, respectively. 



Isokoski et al.: m09bnl0 



Table 3: Source radii and angular diameters for r=100 K. 



Sources 



IRAS20126+4104 

IRAS18089-1732 

G31.41+0.31 



-^T=100K 

[AU] 



tfS,T=100K 
["] 



1753 
2750 
7840 



2.2 
2.4 
2.0 



where the summation is over each line i of the considered 
species. v' is the rest frequency of the line and (/>' is the profile 
function of the line, given by: 



1 



AvV2tt 



-(v-vtf /IAv 2 



(8) 



where Av is the line width in frequency units at 1 /<? height. Av 
can be expressed as a function of the line FWHM in velocity 
units AV by: 



Av: 



V8ln2 



-AV. 



(9) 



The input parameters in the model are the column density 
Ns in cirT 2 , excitation temperature T ex in K, the source diame- 
ter 9s in arcseconds, offset velocity from the source LSR (Local 
Standard of Rest) in km s _1 and the line FWHM in km s _1 . All 
parameters excluding the source diameter are optimized man- 
ually to obtain the best agreement with the observed spectra. 
The source diameter for emission from cold species is allowed to 
vary over a larger area; generally 14" is used. The emission from 
hot-core molecules targeted in our work is assumed to originate 
from the central region with Tjust >100 K. The source diameter 
for the warm emission is calculated using a relation derived from 
dust modeling of a large range of sources (BIS07): 



R 



T=100 K 



2.3 x 10 14 




(10) 



where L/L Q is the luminosity of the source relative to the solar 
luminosity. Table 3 gives the calculated /?t=ioo k radii and diam- 
eters for the observed sources. 

In the analysis for individual molecules, the initial values for 
T ex were based on the T mt from the RTD analysis in case of 
optically thin species. For optically thick species the value of 
T ex from the 13 C-isotope was used. If T ex for the isotopologue 
could not be obtained, the initial temperature was guessed. The 
r ex value was then optimized visually based on the relative in- 
tensity of the emission lines. The simulated emission was not 
allowed to exceed the emission of optically thin lines in the 
spectrum in any of the observed spectral ranges. Coinciding and 
blended transitions, which together contribute to an optically 
thick line, are excluded in the analysis. 

Due to the visual optimization and the possibility of overlap- 
ping lines (particularly in the line-rich source G31), the resulting 
r ex values are only a rough estimate (±50-100 K) and do not dif- 
fer significantly from those from the RTD method. The column 
density for each species was constrained by optically thin, un- 
blended lines, where available. 

For the specific case of CH3OH, which requires two tem- 
perature components for a proper fit, the analysis was also done 
with the CASSIS analysis package 8 . CASSIS has the advantage 



8 CASSIS has been developed by IRAP-UPS/CNRS (http:// 
cassis . irap . omp . eu). 



that it can properly model the emission from overlapping opti- 
cally thick lines, as well as from nested regions of emission. For 
CH3OH, warm emission from the compact inner region may be 
absorbed by the surrounding colder gas, which can influence the 
derived hot core abundances. 



4. Results 

4.1. General results and comparison between sources 

The observed sources, IRAS20126, IRAS18089 and G31, dif- 
fer from each other significantly in the observed chemical rich- 
ness in the JCMT single-dish data. Fig. 2 presents two frequency 
ranges with lines from several observed species. For G31, strong 
lines of all complex organic species are detected, whereas for 
IRAS20126 many targeted lines are below the detection limit. 
Many complex molecules are also found in IRAS 18089 but with 
weaker lines than for G3 1 . Among the serendipitous discoveries 
acetone, CH3COCH3, is possibly seen in G31 and IRAS 18089 
for the first time (see §4.5). Integrated line intensities for all de- 
tected lines and selected upper limits are given in Tables A.l- 
A.14. The rotation diagrams for the detected species are shown 
in Figures B.1-B.13 whereas the optimized parameters in the 
Weeds model for each molecule and source are available online 
in Tables C.l -C.3. 



4.2. Optical depth determinations 

To assess the importance of optical depths effects, the ratio 
of lines of different isotopologues are compared. The expected 
12 C/ 13 C isotope ratio depends on the galactocentric radius, Rgc, 
according to Eq. 1 1 (Wilson & Rood, 1994). 



2 C/ 13 C= (7.5 ± 1.9)R CC [kpc] + (7.6 ± 12.9) 



(11) 



The isotope ratios derived for our sources using Eq. 1 1 are given 
in Table 1 . The galactocentric radii were calculated trigonomet- 
rically from the galactic coordinates, using the IAU value for the 
distance to the Galactic center Rq = 8.5 kpc (Kerr & Lynden- 
Bell 1986). 

Table 4 lists the observed isotopologue intensity ratios for the 
most abundant species in our sources. The CH30H/ 13 CH30H ra- 
tios are derived from a low-energy transition 22,o,+o - 3i,3,+o with 
£ up =44.6 K. No high E up transitions were reliably detected 
for 13 CH30H in our sources in our standard settings and the 
low E up ratios are therefore taken to apply to cold methanol. 
For G3 1.41 +0.31, one additional observation was carried out 
to cover a transition with £ up >100 K. A CH 3 OH/ 13 CH 3 OH in- 
tensity ratio of 4.0 is derived for a transition with E up x2lQ K 
(13i42,-o - 13o,i3,+o)- This indicates that also warm methanol is 
optically thick in this source. 

The CH 3 CN/CH 3 13 CN ratios are derived from the 13 3 - 
12 3 line intensities for G31 and IRAS 18089 and indicate that 
CH3CN is optically thick in these two sources, but not in 
IRAS20126. For H 2 CO and HNCO, isotopologue lines are de- 
tected but for different transitions than the main isotopologues. 
Thus, a model is needed to infer the optical depths. Fits to each 
of the isotopologues independently at a fixed temperature of 
150 K using the /?t=ioo k source size gives column density ra- 
tios that are significantly smaller than the overall isotope ratios, 
suggesting that these species are also optically thick for G31 and 
IRAS 18089. In practice, we have excluded the optically thick 
lines (as indicated by the Weeds model) from the RTD fits for all 
species. 
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Fig. 2: Spectral ranges 238.83-239.26 GHz and 339.94-340.18 GHz with lines from several targeted species for the three sources. 

Table 4: Isotopologue line intensity ratios in the observed sources. Lower limits are those for which the I3 C-isotopologue was not 
detected. 



Species 


IRAS20126+4104 


IRAS18089-1732 


G31.41+0.31 


CH 3 OH/ lj CH 3 OH 
CH 3 CN/CH 3 13 CN 


>6 
>65 


>18 

4.7 


6 
4.5 



4.3. Temperatures 

Table 5 summarizes the derived temperatures from the RTD fit 
for the various species. As also found by BIS07, molecules can 
be classified into cold (<100 K) and warm (>100 K), and our cat- 
egorization is similar to theirs. The Weeds analysis is consistent 
with this grouping. There is however variation in temperatures 
within the groups, with warm species having rotation tempera- 
tures from 70 to 300 K, and cold species from 40 to 100 K. Some 
variation is seen in rotation temperatures of individual species 
between different sources; the rotation temperatures are gener- 
ally higher in G31 than in IRAS18089, while IRAS20126 has 
the lowest of the three. 

The r rot value for CH 3 OH is -300 K for G31 and 
IRAS18089 and -100 K for IRAS20126. For G31, several lines 
with £u P >400 K are detected, which makes the RTD fit more ro- 
bust. For IRAS 18089 and IRAS20126 the accuracy of the RTD 
fits suffers from the small range of E ap in the detected transitions. 
In addition to optically thick lines, low-£up lines have been ex- 
cluded from the fit. These lines are underestimated by the RTD 
fit and probably originate from a colder extended region also 
seen in the 13 C lines (Fig. 3). The T mt from the RTD analysis 
therefore represents the warm CH3OH alone. See Sect. 4.4.1 for 
a more detailed discussion on the CH3OH emission. 

For CH3CN, the T ml values range from ~200 and ~350 K for 
the three sources. The CH3 13 CN rotation diagram gives a value 
of r iot of only 53 K, however. (Fig. 3). This illustrates the large 
uncertainties at high temperatures for optically thick species and 
the possibility of a cold component in addition to the warm one. 

Contrary to the general trend, the r rot value for H2CO is 
somewhat higher (204 K) in IRAS18089 than in G31 (157 K). 
The discrepancy could be influenced by the small number of 
lines used for the analysis. All fitted lines are those belonging 
to the para-H2CO, and the T rot fits are thus not affected by fluxes 
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Fig. 3: Rotation diagrams for 13 CH30H (left panel) and 
CH3 13 CN (right panel) in G3 1 . Triangles represent blended lines 
and are not included from the fit. 



from different spin states. A single transition of H2 13 CO (3^2 - 
2 1,1) was covered, and no information on the excitation temper- 
ature can therefore be obtained from the minor isotopologue. 

The T mt values of the other complex species, HNCO, 
C2H5OH, C2H5CN, NH 2 CHO and CH3OCH3 are around 100 K 
and are slightly higher in G31 than in IRAS 18089. The r rot of 
HCOOCH3 stands out in both sources, in G31 being closer to 
200 K. No lines belonging to these species were detected in 
IRAS20126. 

The species classified as cold by BIS07, CH2CO and 
CH3CCH, indeed show cold rotation temperatures in all sources. 
Not enough lines of CH3CHO or HCOOH, which are also clas- 
sified as cold in BIS07, are observed in our sources for making 
rotation diagrams. 

4.4. Column densities 

Table 7 presents the column densities derived from the 
RTD analysis, Weeds or CASSIS model, and from the 13 C- 
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Table 5: Temperatures (K) derived from the RTD analysis and the Weeds or CASSIS (CH3OH) model. The species are classified 
warm and cold according to BIS07. Square bracketed values are T rot values for optically thick species and round bracketed values 
are T ex values assumed based on temperatures derived from the other sources in this study. 





IRAS20126+4104 


IRAS 18089- 1732 


G31.41+0.31 




RTD 


model 


RTD 


model 


RTD 


model 








warm 


species 






H 2 CO 


123+21 


150 


204±82 


(150) 


157+44 


(150) 


CIROH 


122+17 


300, 14±1 


291±37 


300, 15+2 


323+34 


200, 14+2 


C 2 H 5 OH 


- 


(100) 


85+18 


150 


120+15 


100 


HNCO 


- 


(200) 


92±25 


200 


111+32 


200 


NH 2 CHO 


- 


300 


72±28 


100 


94+50 


300 


CIRCN 


217±352 


(200) 


[346+106] 


(200) 


[311+68] 


(300) 


C 2 IRCN 


- 


(80) 


84+33 


80 


105+23 


80 


HCOOCH, 


- 


(200) 


118+20 


200 


174+11 


300 


CH 3 OCH 3 


- 


(100) 


66+11 


100 


90+6 


100 








cold 


species 






CH 2 CO 


- 


(50) 


71+11 


50 


97+101 


50 


CH3CHO 


- 


(50) 


- 


(50) 


- 


50 


HCOOH 


- 


(40) 


- 


(40) 


- 


40 


CIRCCH 


40±10 


35 


46+12 


40 


67+14 


80 



- means not enough lines were detected for a rotation diagram. 
Typical uncertainties in the Weeds excitation temperatures are ±50 K. 



isotopologues for the optically thick species. Following BIS07, 
the column densities for warm molecules are given as source- 
averaged values (see Eq. 5). The emission from cold molecules 
extends over a larger volume and the values are given as beam- 
averaged column densities. Typical uncertainty of the column 
densities derived from the RTD analysis is ~40 %. 

4.4.1. CH3OH 

An accurate determination of the CH3OH column density is es- 
sential for comparing the abundance ratios of complex organic 
species. For hot-core molecules, it is particularly important to 
quantify the warm CH3OH emission. The column densities of 
CH3OH in BIS07 were determined by the RTD method exclud- 
ing the optically thick lines. The same is done in our analysis. 
Our rotation diagrams however also show emission from low- 
£u P transitions, which are strongly underestimated by the RTD 
fit on the warm lines, providing further evidence for the presence 
of a colder component. We have therefore also excluded these 
transitions. The fit to the higher E up , optically thin lines should 
give the warm CH3OH column density obtained in a similar way 
as BIS07. 

CH3OH emission was also simulated using a two-component 
CASSIS model. A single-component model is not able to simul- 
taneously reproduce the warm and cold lines without overesti- 
mating the lines from intermediate energy levels. Indeed, a better 
agreement is obtained using a two-component model, consisting 
of a warm compact component and a cold extended component. 
Table 6 shows the best model parameters obtained from the fit- 
ting. The warm components are fixed to £>t=iook while the cold 
component is allowed to extend beyond the beam diameter. The 
warm CH3OH column densities derived from the CASSIS fit are 
in agreement with the values derived from the RTD analysis. The 
best fits plotted onto the CH3OH 1 K -6 K transitions (338.5 GHz) 
are shown in Fig. 4. 

Several I3 CH30H lines are detected in G31. However, only 
low-Eup lines are reliably detected since the high-Eup lines are 
very weak or blended. Assigning these lines to the cold com- 
ponent and assuming 7^=20 K (Oberg et al. 2011; Requena- 



- Observed spectrum 
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338.4 338.5 338.6 

Rest Frequency (GHz) 



Fig. 4: CASSIS two-component models for CH3OH emission in 
the 338 GHz region covering the Ik - 6k transitions with K = 
0-6. 



Torres et al. 2008), the beam-averaged 13 CH30H column density 
is 2.1 x 10 15 cm 2 , corresponding to a CH3OH column density 
of 1.7 X 10 17 crrT 2 for the 14" volume. This is more than an or- 
der of magnitude higher than found in the Weeds model for the 
cold component, supporting the optically thick interpretation of 
the cold component. Similarly, for the case of IRAS 18089, the 
detected 13 CH30H lines may indicate a high optical depth. 

As mentioned before, the 13 CH30H line with E up =211 K, 
covered in additional observations for G3 1, reveals optical thick- 
ness in the warm component as well. The CH3OH column den- 
sity derived from this line assuming temperature of 300 K is 2.7 



x 10 18 



cmT 1 and the CASSIS fit as well as the RTD analysis 



underestimate the warm column density by a factor of ~2.5. In 
principle, the CASSIS analysis could be made consistent with 
the 13 CH30H value by letting the warm source size vary as 
well. However, since such information is not known for other 
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Table 6: CASSIS model parameters for CH3OH in the observed sources. 



Source 


Component" 


N s r ex FWHM 9 S Klsr x l 
[cm" 2 ] [K] [kms" 1 ] ["] [kms" 1 ] 


IRAS20126+4104 


comp. 

ext. 


(1.1±0.1)E+17 300 r 8.15±0.60 2.2" -4.05 ± 0.70 
(2.2+^)E+15 14±1 7.00;$*° »14.0 -4.25!^° 


IRAS 18089- 1732 


comp. 

ext. 


(2.0 + ^)E+17 300 r 5A5 +U f 2.4" 33.4±0.2 
(2.4+^)E+15 15±2 3.50+°|° »14.0 32.3+^ 


G3 1.4 1+0.31 


comp. 

ext. 


(1.0±0.2)E+18 200 c 6.35+0.35 2.0" 97.2+0.2 
(1.2±0.6)E+16 14+2 3.95+0.40 »14.0 97.2 + °- 2 



8 "comp." for warm, compact component ; "ext." for colder, more extended component. 

6 Fixed to s , t =iook- 
c Fixed. 



molecules, as well as for consistency with BIS07, we have cho- 
sen to keep the warm source size fixed at the 100 K radius. 

In summary, it seems plausible that at least in our sources, 
methanol emission arises from two temperature components, a 
warm (r ex « 300 K) compact component and a cold (r ex « 20 K) 
significantly more extended component. At least for G31, the 
CH3OH emission is optically thick in both components. 

4.4.2. Other molecules 

Overall, the column densities of the various species follow the 
same pattern in all sources, and hence seem to be well corre- 
lated relative to each other. CH3OH is the most abundant com- 
plex molecule in all sources. The other species have in general 
half to one order of magnitude lower column densities. 

The CH3CN emission is optically thick for G31 and 
IRAS 18089. Column densities from the RTD analysis are there- 
fore underestimated. Due to a lack of sufficient optically thin 
lines, also the Weeds analysis underestimates the column den- 
sities. Values derived from the 13 C-isotopologues are thus more 
reliable, even though some fraction may arise from a colder com- 
ponent. Indeed, the column density for CH3CN derived from the 
13-isotopologue is an order of magnitude larger than that from 
the RTD of the main isotopologue alone. The same procedure 
was used by BIS07. 

For H2CO, the best estimates come from the RTD analysis 
on the optically thin para-H2CO lines. The derived H2CO col- 
umn densities assuming T rot from Table 5 and a statistical ortho- 
to-para ratio of 3 are given in Table 7. The column densities de- 
rived from the I3 C-isotopologue of ortho-P^CO are still larger 
than from RTD analysis using optically thin lines, which indi- 
cates either larger ortho-to-para ratio in these sources (ratio of ~3 
to 5 has been predicted for cold clouds by Kahane et al. 1984) or 
non LTE excitation. BIS07 derived the H2CO column densities 
from the 13 C-isotopologue. 

The column density derived for HCOOCH3 from the RTD 
diagram is significantly smaller than that derived from the Weeds 
model. The HCOOCH3 emission is stronger in the 230 GHz 
beam (20-21") than in the 345 GHz beam (14") probably in- 
dicating significant extended emission. The RTD analysis was 
performed on the entire dataset, while only lines in the 230 GHz 
band was used in the Weeds modeling. 
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Fig. 6: Beam-averaged column densities for the targeted cold 
species. Column densities from the RTD analysis are marked 
in black bars while the red bars show column densities from the 
Weeds model. Arrows indicate upper limits. 
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Fig. 5: Source-averaged column densities for the targeted warm species in &r=ioo k volume. Column densities from the RTD analysis 
are marked in black bars. The red and blue bars show column densities from the Weeds or CASSIS (CH3OH) models and from 
13 C-isotopologue, respectively. Arrows indicate upper limits. 

Table 7: Column densities for the targeted species from the RTD analysis, Weeds or CASSIS (CH3OH) models and those derived 
from 13 C-isotopologues. Column densities for the warm molecules are source-averaged, while those for cold molecules are beam- 
averaged. Column densities where optical depth has a significant effect are labeled as lower limits. 







RAS20126+4104 


IRAS 18089- 1732 




G31.41+0.31 




Species 


RTD 


Model 


N( u C) a 


RTD 


Model 


N( [I C) a 


RTD 


Model 


N( ll cy 


H,CO" 


8.4E+15 


6.0E+15 


<4.7E+16 


1.9E+16 


1.5E+16 


6.3E+16 


6.5E+16 


6.0E+16 


9.6E+16 


CH,OH 


2.9E+16 


1.1E+17 


<4.0E+17 r 


2.2E+17 


2.0E+17 


2.0E+18 C 


1.2E+18 


1.0E+18 


2.7E+18'' 


C 2 H 5 OH 


- 


<0.3E+16 


- 


3.9E+16 


2.5E+16 


- 


1.5E+17 


1.0E+17 


- 


HNCO 


- 


1.0E+15 


<7.0E+15 


2.5E+15 


0.4E+16 


<3.8E+16 


3.6E+15 


0.6E+16 


<8.2E+16 


NH^CHO 


- 


0.3E+15 


- 


8.4E+14 


0.5E+15 


- 


2.1E+15 


0.4E+16 


- 


CH,CN 


1.5E+15 


1.3E+15 


<1.4E+15 


>4.7E+15 


>3.5E+15 


2.2E+16 


>1.6E+16 


>2.0E+16 


1.6E+17 


C 2 H 5 CN 


- 


<0.5E+15 


- 


3.1E+15 


0.4E+16 


- 


2.0E+16 


2.0E+16 


- 


HCOOCH^ 


- 


<0.2E+16 


- 


2.8E+16 


3.0E+16 


- 


1.7E+17* 


1.0E+18 


- 


CH 3 OCH 3 


- 


<0.1E+17 


- 


1.3E+17 


1.0E+17 


- 


5.4E+17 


5.0E+17 


- 


CH 2 CO 


- 


<2.5E+13 


- 


7.9E+13 


7.5E+13 


- 


- 


3.3E+14 


- 


CH3CHO 


- 


<0.5E+14 


- 


- 


<0.5E+14 


- 


- 


5.0E+14 


- 


HCOOH 


- 


<2.5E+13 


- 


- 


<2.5E+13 


- 


3.2E+15 


1.0E+15 


- 


CH,CCH 


2.7E+14 


3.5E+14 


- 


5.2E+14 


0.9E+15 


- 


1.0E+15 


1.4E+15 


- 



* Rotation diagram fit on HCOOCH3 on all lines, while Weeds model on the 230 GHz lines only. 

" N( n C) obtained from the N( n C) adopting a 12 C/ 13 C ratio equal to 70, 54 and 41 for IRAS20126, IRAS18089 and G31.41, respectively. 
Column density from 13 C-isotopologue from ortho-H2 13 CO, corrected using the statistical ortho to para ratio of 3:1. 



c Calculated from the line intensity of transition 2 2 _o.+o 
d Calculated from the line intensity of transition 13 1,1 2,-0 



3i, 3i+ o at 345.084 GHz with £ up =44.6 K, assuming a T ml of 300 K. 



13o 



, at 341.132 GHz with £ up =206 K, assuming a T M of 300 K. 



5. Discussion 

5. 1 . Comparison to massive YSOs without a disk structure 

In order to see the effect of a flattened circumstellar structure 
on the chemistry around the YSO, we compare the temperatures 
and abundances of the complex organic molecules in sources 
studied in this work and those studied in BIS07. In addition to 
our sources (IRAS20126, IRAS18089 and G31), AFGL 2591, 
NGC 7538 IRS1 and G24.78 from the BIS07 sample are now 
known to have disk structures as well (van der Tak et al. 2006; 
Wang et al. 2012; Pestalozzi et al. 2004, 2009; Knez et al. 2009; 
Beltran et al. 2005). A disk-like structure may be present also 
for other sources but current evidence is not as strong as for 
the above sources. For the comparison with the BIS07 sample, 
we use the results primarily from the RTD method to exclude 
method based differences. 



5.1.1. Rotation temperatures 

Figure 7 shows the rotation temperatures of the complex species 
for the observed sources together with those from BIS07. As in 
BIS07, we find that the complex molecules can be divided into 
warm and cold species based on their rotation temperatures. The 
rotation temperatures obtained for the molecules in our sources 
agree with the division; CH 3 OH, H 2 CO, C 2 H 5 OH, HCOOCH3, 
CH3OCH3, HNCO, NH 2 CHO, CH3CN and C 2 H 5 CN show ro- 
tation temperatures generally >100 K while CH 2 CO, CH3CHO, 
HCOOH and CH3CCH have rotation temperatures <100 K. 

The rotation temperatures of the cold molecules show a very 
small scatter (±25 K) from source to source, with or without a 
disk structure, ignoring the one outlier for HCOOH. Of the warm 
species, C 2 H50H has a consistent rotation temperature of 100- 
150 K from source to source. The RTD analysis for these species 
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Fig. 7: Rotation temperatures for selected species in massive YSOs with (open stars) and without (asterisks) observed disk structure. 



is reliable due to low optical depths and lack of anomalous exci- 
tation. 

The other warm species show more scatter in the de- 
rived rotation temperatures. This is particularly the case for 
CH3OH and CH3CN, which have rotation temperatures rang- 
ing from 100 to 350 K. No systematic difference is seen be- 
tween the two source types. Moreover, the HCOOCH3 rotation 
temperatures for AFGL 2591 and NGC 7538 IRS1 are simi- 
lar to the diskless sources, disproving any difference between 
the two source types. The scatter in the rotation temperatures 
could indicate that they exist in environments with different tem- 
peratures, but it may also be caused by optical depth effects 
(H 2 CO, HNCO, CH3OH and CH3CN) and (anomalous) radia- 
tive excitation. The latter has been previously seen for, e.g., 
NH 2 CHO and HCOOCH3 (Nummelin et al. 2000) as well as 
HNCO (Churchwell et al. 1986). For species with significant ra- 
diative excitation, e.g., HCOOCH3, the scatter of the data points 
in the RTD plots is higher and results in additional uncertainty 
since the linear fits do not properly capture this scatter. 



5.1.2. Column densities 

Figures 8 and 9 show the column densities for the targeted 
species in our sources and those from BIS07 for warm and cold 
molecules, respectively. For the warm species the column den- 
sities within the sources vary by 1-2 orders of magnitude. There 
is also a significant variation between sources. G31 is chemi- 
cally richest of the sources, with highest column densities of all 
the targeted species, compared to other sources. IRAS20126 is 
among the chemically poorest sources. The pattern of column 
densities is remarkably similar, however, and the disk sources 
do not stand out. 



5.1.3. Abundance ratios 

Table 8 gives the abundance ratios for the targeted species with 
respect to CH3OH for oxygen-bearing species and HNCO for 
nitrogen-bearing species. Column densities from the RTD anal- 
ysis are primarily used. In cases where column densities are not 
available from the RTD analysis (mainly for IRAS20126) the 
values are taken from the Weeds analysis. 



The resulting abundance ratios are presented in Fig. 10. 
Values from the RTD method are shown in black bars and can 
thus be directly compared to BIS07, whereas those from the 
Weeds model and I3 C-isotopologue are in red and blue bars. 
In general, the abundance ratios in different sources have larger 
variations within the source types than between them. For ex- 
ample, H2CO has both lowest and highest abundance ratio for 
two disk candidates, G31 and IRAS20126, respectively. The 
H2CO results depend somewhat on the analysis method used. 

C 2 H 5 OH, HCOOCH3 and CH3OCH3 abundances peak for 
G31 and IRAS 18089, and are generally lower for other sources. 
CH3OCH3 also peaks for some diskless sources. The N-bearing 
species NH2CHO has a large variation in the abundance ratio 
with respect to HNCO. Among our sources, NH2CHO peaks for 
G31, the source with largest abundance of O-bearing species, 
and the least clear disk structure. 

Figure 1 1 shows the abundance ratios of N-bearing species 
with respect to CH3OH. Again, a lot of variation is seen within 
each source type, with no specific trend found for sources with a 
disk-like structure. 

We stress that the absolute inferred column densities and 
abundance ratios are uncertain by a factor of a few up to an or- 
der of magnitude, as already indicated by the different analysis 
techniques. An independent assessment of the accuracy of the re- 
sults can be made by comparison with the inferred column den- 
sities and abundances of Zernickel et al. (2012) for NGC 6334 I 
(=NGC 6334 IRS), who observed a completely different set of 
lines of the same molecules with Herschel-HIFl and the SMA 
in beams ranging from 2-AO" . Comparison with the results of 
BIS07 shows good agreement within a factor of a few for the 
abundances of several species (CH 3 OH, CH3OCH3, CH3OCHO, 
C2H5OH, C2H5CN) whereas others differ by an order of mag- 
nitude (H2CO, CH3CN). The species that show the largest dis- 
crepancy are those with highly optically thick lines and without 
a large set of isotopologue lines. Thus, a combination of dif- 
ferences in adopted source sizes and optical depth effects can 
account for the discrepancies. Because our approach is the same 
for all sources, the relative values from source to source should 
be more reliable. 
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Fig. 8: Source-averaged column densities for warm molecules. Column densities from the RTD analysis are marked in black bars. 
Sources without disk structure are marked with an asterisk. The red and blue bars show column densities from the Weeds or CASSIS 
(CH3OH) models and from 13 C-isotopologue, respectively. Upper limits are indicated with arrows. 
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Fig. 9: Beam-averaged column densities for cold molecules in sources observed in this study and those from BIS07. Sources without 
disk structure are marked with an asterisk. Column densities from the RTD analysis are marked in black bars while the red bars 
show column densities from the Weeds model. Upper limits are indicated with arrows. 



5.2. Chemical and physical implications 

The main goal of this study is to investigate similarities and 
differences between sources with and without a disk-like struc- 
tures. The presence of a flattened accretion disk should allow 
the UV radiation from the central object to escape more read- 
ily and then impinge on the gas and dust in the outer walls of 
a flared disk or outflow cavity (Bruderer et al. 2010). Increased 
UV radiation could manifest itself as enhancement in the com- 
plex organic species produced through UV photodissociation 



from CH3OH in the solid state, such as CH3OCH3, C 2 H 5 OH 
and CH3OCHO (Oberg et al. 2009). In this scenario, CH3OH is 
photodissociated into various radicals such as CH3, CH3O and 
CH2OH which can become mobile at higher ice temperatures 
(20-40 K) and form the observed complex organic molecules 
(Garrod & Herbst 2006). Higher temperatures favor diffusion 
of larger radicals resulting in the formation of larger complex 
organic molecules compared with small molecules like H2CO. 
Another related parameter that plays a role is the CO content of 
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Table 8: Abundance ratios of complex species in the observed sources and those in other chemically rich environments. The abun- 
dances are column densities with respect to CH3OH for oxygen-bearing species and with respect to HNCO for nitrogen-bearing 
species. The CH3OH and HNCO column densities are taken from the RTD analysis when available. 



Source 



H 2 CO C2H5OH NH 2 CHO HCOOCH3 CH 3 OCH, 



Low-mass outflows 



Low-mass protostars 



Low-mass hot corinos 



Refs. 



Massive protostellar objec 


ts with a 


disk structure 










IRAS20126+4104 


0.29" 


<0.10 


<0.30 


<0.34 


<0.34 


this work 


IRAS 18089- 1732 


0.09" 


0.18 


0.34 


0.13 


0.60 


this work 


G31.41+0.31 


0.06" 


0.12 


0.59 


0.14 


0.46 


this work 


AFGL2591 


0.28" 


<0.02 


<1.00 


0.51 


<0.16 


1 


NGC 7538 IRS 1 


0.21° 


0.05 


0.25 


0.12 


<0.13 


1 


G24.78 


0.23° 


0.03 


0.14 


0.11 


0.43 


1 



Massive protostellar objec 


ts with no 


detected disk structure 








G75.78 


0.20° 


<0.02 


0.05 


0.06 


0.21 


1 


W33A 


0.27° 


0.02 


0.32 


0.13 


0.14 


1 


NGC 6334 IRS 1 


0.13° 


0.02 


0.17 


0.12 


0.60 


1 


W 3(H20) 


0.18° 


0.01 


0.27 


0.05 


0.15 


1 



Bl-b core 


1.2 


<1.0 


2.3 


<0.8 


2 


SMM4-W 


0.6 


<0.4 


3.5 


1.1 


3 


LI 157 outflow 


- 


0.7 


1.8 


- 


4,5 



SMM1 


6.6 


<3.4 


10 


5.3 


3 


SMM4 


2.2 


<0.6 


<1.0 


0.8 


3 


NGC 1333 IRAS 4Aenv. 


- 


- 


56 


<22 


6,7,8 


NGC 1333 IRAS 4B env. 


- 


- - 


26 


<19 


6,7,8 


IRAS 16293 env. 


4 


- 


30 


20 


8, 9, 10 



IRAS 16293 A 


<0.02 


1.4 


0.6 


0.6 


11, 12 


IRAS 16293 B 


0.6 


- 


0.8 


1.6 


11, 12 


NGC1333IRAS2A 


- 


- 


- 


2 


13, 14 



" para-ILCO. 

ortho-H 2 CO. 

References. (1) Bisschop et al. (2007); (2) Oberg et al. (2010); (3) Oberg et al. (2011); (4) Bachiller & Perez Gutierrez (1997); (5) Arce et al. 

(2008); (6) Maret et al. (2005); (7) Bottinelli et al. (2007); (8) Herbst & van Dishoeck (2009); (9) van Dishoeck et al. (1995); (10) Cazaux et al. 

(2003); (11) Kuan et al. (2004); (12) Bisschop et al. (2008); (13) Huang et al. (2005); (14) J0rgensen et al. (2005). 



ices, with H2CO, CHO- and COOH-containing species enriched 
in cold, CO-rich CH3OH ices in which CO has not evaporated. 
Indeed Oberg et al. (2011) show that the increased abundance 
of HCOOCH3 in low-mass YSOs compared with the high-mass 
sample of BIS07 may be due to this effect. Although only lim- 
ited laboratory experiments exist on N-containing molecules, the 
abundance of N-bearing complex organic molecules could be en- 
hanced due to photodissociation of N2 or NH3. 

Our results show that there is no consistent enhancement of 
complex organic species produced through CH3OH photochem- 
istry, nor of N-bearing species, that can be attributed to the pres- 
ence of a flattened disk structure on the scales probed by our 
data. The variations within source types are larger than between 
the source types. One explanation for this lack of differenta- 
tion could be that all high-mass sources have a source structure 
through which UV radiation can impact the larger surroundings. 
In particular, just like low-mass sources, all high-mass YSOs are 
expected to have outflows and cavities through which UV radi- 
ation can escape and affect the chemistry, whether or not they 
have a large disk-like structure. This scenario can be tested with 
future high angular resolution data on <1" scale with ALMA 
which should then reveal the emission from complex molecules 
coating the walls of the outflow cavities. 

An alternative explanation is that enhanced temperature in 
the YSO environment is not needed for the production of com- 
plex organic molecules and that they are already formed in the 
prestellar stage through UV- and cosmic -ray processing of cold 
ices, using just the internal UV field produced by interaction of 



cosmic rays with H2 (Prasad & Tarafdar 1983) rather than that of 
the star. This scenario has gained support over the last few years 
with the detection of HCOOCH3 and CH3OCH3 in cold low- 
mass cores away from the YSO where the molecules are either 
released by shocks (Arce et al. 2008; Oberg et al. 2010) or by 
photodesorption (Bacmann et al. 2012) although the precise for- 
mation mechanisms are still unclear. Gas-phase processes may 
also contribute in some cases. Again, high spatial resolution and 
sensitivity such as ALMA will be key to testing this scenario. 
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Fig. 10: Abundance ratios of complex species with respect to CH3OH for oxygen-bearing species and HNCO for nitrogen-bearing 
species. The black bars indicate abundance ratios calculated with the CH3OH and HNCO column density derived from the RTD 
analysis similar to BIS07. The red and blue bars indicate abundance ratios where CH3OH and HNCO column densities are derived 
from CASSIS or Weeds (HNCO) model, respectively. Upper limits are marked with arrows. 
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Fig. 11: Abundance ratios of N-bearing species HNCO, CH 3 CN, C 2 H 5 CN and NH 2 CHO with respect to CH 3 OH (from RTD anal- 
ysis). Arrows indicate upper limits. 



6. Summary and conclusions 

We have carried out a partial submillimeter line survey, tar- 
geting several complex organic species toward selected mas- 
sive YSOs with strong evidence of circumstellar accretion disks, 
IRAS20126+4104, IRAS 18089-1732 and G3 1.41 +0.31. This is 
the first time that molecular abundances are reported for these 
well-known sources. The analysis is performed using both the 
rotation diagram method and spectral modeling. The inferred 
rotation temperatures and molecular abundances are compared 
to sources without reported disk structures analysed using the 
same techniques. The molecules can be divided into two differ- 
ent groups based on their rotation temperatures, independent of 
source type. In particular, the cold (<100 K) species have re- 
markably constant rotation temperatures from source to source. 



The warm (>100 K) species exhibit more scatter possibly due 
to optical depth effects, non-LTE conditions and radiative ex- 
citation. The column densities peak for the same sources, with 
G3 1 .41+0.3 1 being chemically the richest of the studied sources. 
The relative abundances of species follow the same pattern, 
and no chemical differentiation could be established between 
the two source types within the (considerable) uncertainties. The 
lack of chemical differentiation between massive YSOs with and 
without observed disk-like structure suggests similarity in the 
physical conditions and level of UV exposure for all sources, 
for example through outflow cavities. Alternatively, the complex 
molecules may already be formed in the cold prestellar stage 
under similar conditions. Ultimately and indirectly, these results 
imply that the mechanism of the formation of massive stars does 
not differ significantly from source to source. Future high angu- 
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lar and high sensitivity observations of optically thin lines will be 
key to distinguishing the different scenarios and derive accurate 
abundance ratios that can be directly compared with chemical 
models. 
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Appendix A: Detected lines per species for all 
sources 

The line assignment and detection is based on Gaussian fitting 
with the following criteria: (i) the fitted line position has to be 
within ± 1 MHz of the catalog frequency, (ii) the FWHM is con- 
sistent with the those given in table 1 and (iii) the peak intensity 
has at least a S IN - 3. The errors on the integrated intensities 
are computed as follows. 

The integrated main beam temperatures, J TMBdV, were ob- 
tained by Gaussian fits to the lines (Eq. A.l). 

f T MB dV = c s 'T {) AV (A.l) 

with 

st _ V7r/ln2 

C ~ 2 

where To is the peak intensity and AV is the FWHM of the line. 
The error, d J Tmb^V, is calculated from Eq. A. 2. 



'I 



2 
TO 



d T MB dV = c s '[AVd(T ) + T d(AV)] (A.2) 

with 

d(T ) = -yjrms 2 + (cal x T ) 2 + crl 

and 

d(AV) = cr AV 

where rms is the root mean square amplitude of the noise in the 
spectral bin 6v, cal is the calibration uncertainty of the telescope, 
and cr's are the statistical errors on To and AV from the Gaussian 
fits. 

The errors on the integrated intensities derived from Eq. A.2 
include all statistical errors from the Gaussian fit. 

d f T MB dV = c s 'AV Vrms 2 + (cal x T ) 2 (A.3) 

For undetected transitions the upper limits were determined 
as 3cr limits (Eq. A. 4) using: 

cr= 1.2 Vc5vA Vrms, (A.4) 

where 1 .2 is the coefficient related to the calibration uncertainty 
of 20%. 

CLASS was used to determine the Gaussian fits and the un- 
certainties in the individual parameters. The formal errors on the 
integrated intensities derived from Eq. A.3 in some cases yield a 
S /N < 2.5. This is caused by: (i) a conservative estimate of the 
statistical error on the FWHM parameter in the Gaussian fitter 
of CLASS; (ii) all statistical errors are included into our error 
calculation. Considering the higher S /N on the peak intensity, a 
more traditional error estimate without the statistical errors from 
the Gaussian fit (Eq. A.3) would result in a S /N > 3 for the in- 
tegrated intensity as well. All weak line fits were confirmed by 
visual inspection. 
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Table A.l: Observed line fluxes J Tmb^V (K km s -1 ) for H2CO and its isotopic species. 



Frequency 


Transition 




A 


Sources 






[GHz] 




[K] 


[s- 1 ] 


IRAS20126+4104 


IRAS18089-1732 


G3 1.4 1+0.31 


H 2 CO 


218.222^ 


3o,3 - 2o,2 


21 


2.8(-4) 


2.82 (0.65) 


- 


>7.14(1.71) 


218.476^ 


3 2,2 - 22,1 


68 


1.6(-4) 


- 


- 


4.15(1.03) 


363.946^ 


52,4 - 42,3 


100 


1.203) 


4.65 (1.08) 


9.92(2.15) 


>20.29(3.31) 


364.103 P 


54,1/2 -44,0/1 


241 


5.0(4) 


1.90(0.58) 


4.25 (1.07) 


10.97 (2.88) 


364.275 ° 


53,3 - 43,2 


143/158 


8.9(4) 


>5.08(1.50) 


>12.62(3.34) 


>21.56(5.08) 


364.289 ° 


5 3 ,2-4 3 ,i 


143/158 


8.9(4) 


>5.78 (1.71) 


>10.68 (2.83) 


>20.46 (4.85) 


H 2 13 CO 


219.909 " 


3i,2 - 2ii 


33 


2.6(4) 


<0.34 


1.08(0.60) 


1.60 (0.55) 



The notation a(-b) stands for a x 10 

> means lower limit (optically thick line) 

- means frequency not observed 

P para-H 2 CO 

ortho-H 2 CO 
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Table A. 2: Observed line fluxes J TmbcIV (K km s _1 ) for CH3OH and its isotopic species. 



The notation a(-b) stands for a x 10~* 

" from additional observations on G3 1 only 

* line close to the edge of the frequency band 
c used to constrain the RTD fit 

> means lower limit (optically thick line) 

< means upper limit 

~ means uncertain detection, S/N < 2 

- means frequency not observed 

: means blended frequency 

:: cold (£ up <100 K) CH3OH line (treated as a blend with cold component) 

* calculated flux (see Eq. 4) 



Frequency 


Transition 


f 


A 


Sources 






[GHz] 




[K] 


[s- 1 ] 


IRAS20126+4104 


IRAS18089-1732 


G31.41+0.31 


CH3OH 


218.440 


4-2,0 _ 3 1.0 


45 


6.9(-5) 


- 


- 


10.10(2.22):: 


219.984 


253.o _ 244.0 


802 


3.0(-5) 


<1.33 


<0.24 


<0.45 (0.24): 


219.994 


235.o - 22 6 .o 


776 


2.6(-5) 


<1.33 


<0.24 


< 1.74 (0.56): 


240.739 


263,-0 - 262.+0 


864 


1.1 (-4) 


<0.20 


<0.40 (0.32) 


<0.78 (0.47): 


240.818 


5l.2-4,, 2 


834 


8.5(-5) 


<0.20 


<0.23 


<3.58 (0.88): 


240.861 


5-4,2 ~ 4_4.2 


779 


3.2(-5) 


<0.20 


<0.27 


<0.81 (0.37): 


240.870 


5o,2 - 4o,2 


769 


8.9(-5) 


<0.20 


<0.90 (0.71): 


<3.26(1.06): 


240.916 


53,-2/+2 - 43 -2/+2 


693 


5.7(-5) 


<0.20 


<1. 17 (0.46): 


1.43 (0.45) 


240.932 


54,-2/+2 - 44 ,-2/+2 


649 


3.2(-5) 


<0.20 


<0.57 (0.20) 


-0.42 (0.21) 


240.937 


5-2,2 "" 4-2,2 


680 


7.4(-5) 


<0.20 


<1.03 (0.69): 


-0.49 (0.24) 


240.939 


5o,+2 ~ 4 ,+2 


543 


1.2(-4) 


<0.20 


<1.03 (0.69): 


2.54 (0.86) 


240.948 


5 3 ,2 - 4 3 ,2 


656 


5.6(-5) 


<0.20 


-0.24(0.12)* 


-0.74 (0.30)* 


240.952 


52,2 - 42,2 


621 


7.5(-5) 


<0.20 


-0.24(0.12)* 


-0.74 (0.30)* 


240.959 


5-1,2-4-1.2 


567 


8.5(-5) 


<0.20 


1.18(0.42): 


1.74(0.61): 


240.961 


5i.+i -4[.+ i 


360 


8.5(-5) 


<0.20 


1.18(0.42): 


1.74(0.61): 


241.043 


22-e.o — 23-5.o 


776 


3.4(-5) 


- 


-0.35 (0.48) 


<0.38 (0.54): 


330.794 


8-3,0 - 9-2,0 


146 


8.0(-5) 


-0.28 (0.57) 


3.23(1.12) 


8.06 (2.43) 


331.220 


16-i,o - 15-2,o 


321 


7.8(-5) 


<0.39 c 


2.66 (0.79) 


>3. 19 (0.91)* 


338.124 


7o,o - 60,0 


78 


2.5(-4) 


- 


- 


30.57 (11.11):: 


338.345 


7-i,o - 6_i.o 


70 


2.5(-4) 


7.29(2.05):: 


10.34 (2.58):: 


38.03 (9.06):: 


338.405 


76,0 - 66,0 


244 


6.7(-5) 


<8.90 (2.54): 


<12.95 (3.14): 


<44.10 (10.87): 


338.409 


7o,+o - 6o,+o 


65 


2.5(-4) 


8.90(2.54):: 


12.95(3.14):: 


44.10(10.87):: 


338.431 


7-6,0 - 6-6,0 


254 


6.7(-5) 


-0.72(0.14) 


-2.23 (1.06) 


8.49 (2.90) 


338.442 


76.-0/+0 - 66 -0/+0 


259 


6.7(-5) 


-1.20(0.79): 


-2.29 (1.27) 


8.80 (3.44) 


338.457 


7-5,0 - 6-5,0 


189 


1.2(-4) 


1.80(0.66): 


-2.55 (1.44) 


10.12(3.74) 


338.475 


7s,o - 65,0 


201 


1.2(-4) 


-1.43 (0.94): 


-1.89(1.21) 


9.22 (3.46) 


338.486 


75,-o/+o - 65 -0/+0 


203 


1.2(-4) 


-2.86 (0.94): 


>3.40 (3.00) 


>13.31 (5.28) 


338.504 


7-4,0 - 6-4,0 


153 


1.7(-4) 


-1.78(0.85) 


-3.39(1.88) 


> 12. 14 (4.99) 


338.5126 


74/4 -0/+0 - 64/4 -0/+0 


145 


1.7(-4) 


>2.61 (0.94): 


>5.87(1.81): 


>26.88 (9.59): 


338.5129 


72,-0 - 62,-0 


103 


2.3(-4) 


>2.61 (0.94): 


>5.87(1.81): 


>26.88 (9.59): 


338.530 


74.0 - 64,0 


161 


1.7(-4) 


-1.34(0.59) 


-3.69 (2.22) 


11.33(4.40) 


338.541 


73.+0/-0 - 63.-1-0/-0 


115 


2.1 (-4) 


>3.75 (1.24) 


>8.71 (2.71) 


>29.00 (8.01) 


338.560 


7-3,0 - 6-3.0 


128 


2.1 (-4) 


-2.14(0.95) 


>4.50 (1.87) 


14.94 (4.70) 


338.583 


73,0 - 63,0 


113 


2.1 (-4) 


2.74(1.17) 


>4,02(1.61) 


18.35(5.14) 


338.615 


7i,o -61,0 


86 


2.5(-4) 


4.48 (1.33):: 


9.01 (2.57):: 


28.49 (6.87):: 


338.640 


72.+0 - 62,+0 


103 


2.3(-4) 


>2.99 (1.07) 


>6.46 (2.36) 


21.12(5.79) 


338.722 


72,0 - 62,0 


87 


2.3(-4) 


2.60(0.70):: 


5.24(1.36):: 


17.99 (8.33):: 


338.723 


7-2,0 _ 6_2,0 


91 


2.3(-4) 


2.60(0.70):: 


5.24(1.36):: 


17.99(8.33):: 


340.141 


22,+o - 3i,+o 


45 


4.1 (-5) 


1.56(0.71):: 


3.31 (1.04):: 


12.35 (2.83):: 


342.730° 


13l, 12,-0 - 13o,13,+0 


147 


2.4(-5) 


- 


- 


14.19(3.54) 


344.971 


127,1-116.1 


762 


1.3(-4) 


<0.12 


- 


-2.27(1.07) 


364.159 


9 3 ,i-9 2 ,i 


522 


6.5(-5) 


-0.20 (0.41) 


2.97(1.10): 


2.66 (0.77) 


1J CH 3 OH 


222.468 


21l,20,0 -21o,21,0 


544 


3.6(-5) 


<0.23 


<0.32 


<0.59 


225.404 


162,14,0 - 153,12,0 


331 


2.0(-5) 


<0.19 


<0.45 (0.35): 


<3.49(1.55): 


338.760 


13o,13,+0 - 12l,12,+0 


206 


2.2(-4) 


<0.20 


<0.39 


-5.55 (2.36) 


341.132° 


13l, 12,-0 - 13o,13,+0 


211 


2.4(-5) 


- 


- 


3.55 (2.22) 


345.084 


22.0.+0 - 3l,3,+0 


45 


2.9(-5) 


<0.29 


<0.24 


<1.69 (0.77): 


345.133 


4o,4,0 - 3-1.3.0 


36 


8.2(-5) 


<0.29 


1.81 (0.70) 


4.06(1.25) 


354.446 


4l,3,0 - 3o,3,0 


44 


1.3(-4) 


<0.18 


1.68(0.65) 


6.32(1.59) 


363.924 


143,12,1 _ 142.13,1 


655 


4.7(-5) 


<0.24 


<2.06 (0.94): 


< 1.40 (0.92): 
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Table A.3: Observed line fluxes j T MB dV (K km s _1 ) for C 2 H 5 OH. 



Frequency 


Transition 


'-'up 


A 


Sources 






[GHz] 






[K] 


[s- 


'] 


IRAS20126+4104 


IRAS18089-1732 


G3 1.4 1+0.31 


C2H5OH 


218.461 


5 3 ,2,2 " Z 


2,3,2 


24 


6.6 


-5) 


- 


- 


<2.60 (0.52): 


218.554 


215,16,2 " 


-214,17,2 


226 


6.2 


-5) 


- 


- 


-0.77 (0.65) 


222.217 


20 5 ,,5,2 " 


~ 2O4.16.2 


208 


6.5 


-5) 


<0.22 


-0.56 (0.59) 


< 1.29 (0.26): 


222.419 


204,17,1 " 


" 195,15,0 


256 


1.1 


-5) 


<0.22 


<0.34 (0.30): 


< 1.1 8 (0.24): 


222.519 


26 4 _23,2 " 


" 263.24,2 


316 


6.7 


-5) 


<0.22 


<0.27 


-0.26 (0.24) 


225.105 


13l0,3/4,C 


- 12l0,2/3,0 


255 


4.2 


-5) 


<0.19 


<0.25 


-0.66 (0.42) 


225.108 


13ll,2/3,C 


- 12i 1,1/2,0 


280 


2.9 


-5) 


<0.19 


<0.25 


-0.40 (0.42) 


225.110 


136,8/7,1 


- 126,7/6,1 


181 


8.0 


-5) 


<0.19 


<0.25 


1.12(0.42) 


225.112 


139,5/4,0 


- 129,4/3,0 


231 


5.3 


-5) 


<0.19 


<0.25 


-0.74 (0.42) 


225.116 


13l2,l/2,C 


- 12l2,0/l,0 


308 


1.5 


-5) 


<0.25 


<0.25 


<0.55 (0.43): 


225.131 


138.5/6,0 


_ 128,4/5,0 


211 


6.4 


-5) 


<0.19 


<0.25 


-1.16(1.12) 


225.171 


137,6/7,0 


- 127,5/6,0 


192 


7.3 


-5) 


<0.19 


<0.25 


-1.86(0.84) 


225.210 


19 5 ,14,2 " 


" 19 4 ,15,2 


191 


6.7 


-5) 


<0.19 


<0.25 


-1.98(1.05) 


225.229 


172,15,2 " 


- 163.14,2 


137 


2.5 


-5) 


<0.19 


-0.37 (0.24) 


<2.76 (0.55): 


225.249 


136.8/7,0 


_ 126,7/8,0 


176 


8.1 


-5) 


<0.19 


-0.40 (0.30) 


-2.34 (0.93) 


225.279 


135,9,1 - 


125,8,1 


168 


8.7 


-5) 


<0.19 


<0.25 


-1.61 (0.77)* 


225.283 


135,8,1 - 


125,7,1 


168 


8.7 


-5) 


<0.19 


<0.25 


-1.61 (0.77)* 


225.400 


135.9.0 - 


125,8,0 


163 


8.7 


-5) 


<0.19 


-0.13(0.15) 


0.74 (0.32)* 


225.404 


135.8,0 - 


125,7,0 


163 


8.7 


-5) 


<0.19 


-0.13(0.15) 


0.74 (0.32)* 


225.457 


133,11,1 " 


- 123,10,1 


148 


9.7 


-5) 


<0.19 


<0.25 


-0.95 (0.48) 


238.841 


212,19,0 - 


-21 1,21,1 


258 


4.8 


-5) 


- 


<0.22 


<2.28 (0.46): 


239.020 


284.25,2 " 


- 275.22,2 


362 


2.9 


-5) 


<0.21 


<0.22 


<4.79 (0.96): 


239.186 


16o.l6,l - 


" 15l,l4,0 


171 


4.4 


-5) 


<0.21 


<0.22 


<0.92(0.18): 


240.654 


42,2,1 - : 


1,2,0 


75 


4.1 


-5) 


- 


<0.22 


- 


240.782 


152,13,0 " 


- 14 U3 ,i 


163 


2.9 


-5) 


<0.2 


<0.23 


<3.38 (0.68): 


240.839 


14l,13,0 - 


- 13o,13,l 


147 


5.1 


-5) 


<0.20 


-0.20 (0.33) 


-0.85 (0.52) 


330.985 


63,4,1--' 


2,4,0 


90 


1.2 


-4) 


<0.31 


- 


-1.46(0.84) 


331.027 


135,8,1 - 


134,10,0 


168 


1.1 


-4) 


<0.29 


<0.29 


<1.84 (0.37): 


331.079 


43 7 ,37,2 " 


-43 6 . 38 .2 


860 


2.2 


-4) 


<0.29 


-0.36 (0.39) 


-0.15 (0.22) 


331.095 


277.20,2 " 


" 276,21,2 


380 


2.1 


-4) 


<0.29 


<0.29 


-0.33 (0.55) 


338.088 


25 1,24,1 - 


" 24 2 ,22,0 


332 


5.9 


-5) 


- 


- 


-3.64(1.99) 


338.099 


187,11,2 - 


- 186,12,2 


205 


2.1 


-4) 


- 


- 


-2.05 (1.54)* 


338.110 


187,12,2 - 


" 186,13,2 


205 


2.1 


-4) 


- 


- 


-2.05 (1.54)* 


338.163 


102,8,1 " 


9l,8,0 


113 


8.6 


-5) 


- 


- 


2.05 (0.98) 


338.412 


177,10,2 " 


" 176,11,2 


190 


2.1 


-4) 


<4.41 (1.31): 


<13.78(3.29): 


<21.70(5.13): 


338.417 


177,11,2 " 


- 176,12,2 


190 


2.1 


-4) 


<4.41 (1.31): 


<0.71 


<21.70(5.13): 


338.672 


167.9,2 - 


166,10,2 


176 


2.0 


-4) 


<0.32 


<0.83 


1.63 (0.85)* 


338.674 


167.10,2 - 


- 166,11,2 


176 


2.0 


-4) 


<0.32 


<0.83 


1.63 (0.85)* 


339.979 


9 4 ,6,2 - i 


3,5,2 


58 


2.2 


-4) 


<0.15 


-0.55 (0.47) 


3.23 (0.93) 


345.174 


77,0/1,0 - 


66,0/1,1 


140 


2.5 


-4) 


<0.26 


-1.16(0.70) 


-2.28(1.56) 


345.229 


21 1,21,0 ■ 


- 20 1.20,0 


242 


3.7 


-4) 


<0.29 


-0.58 (0.37) 


<5.52 (1.10): 


345.295 


21 1,21,1 - 


- 20 1.20,1 


246 


3.7 


-4) 


<0.29 


-1.03 (1.00) 


-2.91 (1.31) 


345.312 


223J9.1 " 


-214,17,0 


286 


5.0 


-5) 


<0.29 


<0.19 


<0.39 


345.333 


21o,2i,o - 


- 20o.20,0 


242 


3.7 


-4) 


<0.29 


<6.54 (1.31): 


<17.35 (3.47): 


345.408 


21o,2i,i - 


- 20o,2o,i 


246 


3.7 


-4) 


<0.29 


-0.75 (0.47) 


-2.47(1.13) 


352.858 


21 1,20,2 " 


- 202,19,2 


196 


1.2 


-4) 


<0.18 


-0.51 (0.43) 


-1.56(0.78) 


353.034 


12^ 9 7 — 


1 12,10,2 


77 


1.9 


-4) 


<0.16 


-0.61 (0.59) 


3.12(1.00) 


354.758 


20 3 ;, 7 ", - 


- 19 3 ,.6,1 


249 


4.0 


-4) 


<0.18 


-0.59 (0.52) 


3.96(1.05) 


363.968 


217,15/14,1 -207,14/13,1 


314 


3.9 


-4) 


<0.22 


-0.64 (0.44) 


2.66(1.02) 


364.001 


218,14/13.0 - 208,13/12,0 


327 


3.7 


-4) 


<0.22 


-0.70 (0.48) 


2.28 (0.80) 


364.233 


217,15/14,0 -207. 14/13.0 


309 


3.9 


-4) 


<0.24 


-0.64 (0.56) 


-1.94(1.48) 



The notation a(-b) stands for a x 10" 
< means upper limit 

- means uncertain detection, S/N <,1 

- means frequency not observed 
: means blended line 

* calculated flux (see Eq. 4) 
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Table A. 4: Observed line fluxes J TMBdV (K km s _1 ) for HNCO and its isotopic species. 



Frequency 
[GHz] 



Transition 



p 

'-'Up 

[K] 



A 

[s- 1 ] 



Sources 
IRAS20126+4104 



IRAS18089-1732 G31.41+0.31 



HNCO 
















219.657 


103,8/7,9/11/10 


_ 93,7/6,9/10/8 


447 


1.4(-4) 


<1.33 


<0.73 (0.70): 


<1.52(0.63): 


219.736 


102,9/8,9/11/10 


_ 92,8/7,9/10/8 


231 


1.4(-4) 


<0.28 


<0.84 (0.75): 


-0.64 (0.91) 


219.798 


10o,10,9/ll/10 - 


9o,9,9/10/8 


58 


1.5(-4) 


<0.28 


1.31 (0.50) 


1.22(0.55) 


240.876 


Hi, 11,10/12/11 


- 10l, 10,10/11/9 


113 


2.004) 


<0.20 


<0.82 (0.54) 


<3.10(0.93): 


330.849 


15l, 14,14/16/15 


- 14l, 13,14/15/13 


170 


5.204) 


-0.32(0.31) 


<1.53 (0.71): 


2.01 (0.71) 


352.898 

TT*t13^/-> 


16l, 15. 15/17/16 


- 15l, 14,15/16/14 


187 


6.304) 


<0.18 


-1.84(1.00) 


<3. 12 (1.02): 



219.664 


103,8/7,10 - 93,7/6,10 


47 


1.4(-4) 


<1.32 


<0.21 


<0.41 


219.740 


102,9,11/9/10 - 92,8,10/8/9 


31 


14(4) 


<0.28 


<0.24 


<0.26 


219.744 


102,8,11/9/10 - 92,7,10/8/9 


31 


14(4) 


<0.28 


<0.24 


<0.26 


219.804 


10o,10,10/9/8 - 9o,9,10/9/8 


58 


1.5(4) 


<0.28 


<0.24 


<0.63 


240.881 


1 111, 11, 12/11/10 - 10l0.ll, 11/10/9 


113 


2.0(-4) 


<0.20 


<0.23 


<0.98 


330.860 


15l, 14,16/14/15 - 14l, 13,15/13/14 


170 


5.2(4) 


<0.31 


<0.32 


<0.31 



The notation a(-b) stands for a x 10 
< means upper limit 

- means uncertain detection, S/N < 2 

- means frequency not observed 
: means blended line 



Table A.5: Observed line fluxes J T MB dV (K km s 1 ) for NH 2 CHO. 



Frequency 


Transition 


F 

'-'Up 


A 


Sources 






[GHz] 




[K] 


[s" 1 ] 


IRAS20126+4104 


IRAS18089-1732 


G31.41+0.31 


NH 2 CHO 


218.460 


10l0,l,9 -99,1,8 


61 


7.5(-4) 


- 


- 


<2.45 (0.49): 


339.716 


16l6,8,8/9 - 15l5,8,7/8 


329 


2.2(-3) 


<0.15 


<0.17 


<4.04(0.81): 


339.781 


16l6,7,10/9 - 15l5,7,9/8 


284 


2.3(-3) 


-0.74 (0.42) 


-0.29 (0.29) 


<2.66 (1.20): 


339.904 


16l6,6,ll/10 ~ 15l5,6,10/9 


246 


2.5(-3) 


<0.89 (0.39): 


-0.48 (0.30) 


-1.41 (0.76) 


340.135 


16l6,5,12 - 15 15,5,1 1 


213 


2.6(-3) 


<0.78(0.16): 


< 1.40 (0.28): 


<1. 11 (0.22): 


340.139 


16l6,5,ll - 15l5,5,10 


213 


2.6(-3) 


<0.78(0.16): 


< 1.40 (0.28): 


<1. 11 (0.22): 


345.183 


17n,0,17 - 16l6,0,16 


151 


3.0(-3) 


<0.20 


1.16(0.53) 


2.42 (0.94) 


345.327 


16l6,1.15 - 15l5,l,14 


145 


3.0(-3) 


<0.43 


< 1.87 (0.37): 


<5.46(1.09): 



The notation a(-b) stands for a x 10 
< means upper limit 

- means uncertain detection, S/N <2 

- means frequency not observed 
: means blended line 
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Table A. 6: Observed line fluxes J T^BdV (K km s -1 ) for CH3CN and its isotopic species. 



Frequency 
[GHz] 



Transition E up 
[K] 



A 
[s- 1 ] 



Sources 
IRAS20126+4104 



IRAS18089-1732 G31.41+0.31 



CH 3 CN 
















238.844 


13 s - 


-12* 


537 


5.0(-4) 


- 


<0.22 


<0.17 


238.913 


13,- 


-12, 


430 


5.8(-4) 


<0.21 


-0.37 (0.27) 


1.82(0.62) 


238.972 


13 6 - 


-12 6 


337 


6.4(-4) 


-0.33 (0.38) 


-0.84 (0.37) 


3.22(1.11) 


239.023 


13 5 - 


-12 5 


259 


6.9(-4) 


-0.93 (0.49) 


<6.43(1.29): 


3.19(0.88) 


239.064 


13 4 - 


-12 4 


195 


7.4(-4) 


-0.51 (0.38) 


1.17(0.45) 


5.25 (1.44) 


239.096 


13 3 - 


-12 3 


145 


7.7(-4) 


1.04(0.40) 


1.64(0.55) 


5.21 (1.37) 


239.120 


13 2 - 


-12 2 


109 


7.9(-4) 


-0.63 (0.38) 


1.80(0.61) 


5.71 (1.37) 


239.133 


13i- 


-12i 


87 


8.1(-4) 


0.73 (0.26)* 


1.72(0.52)* 


4.56(1.19) 


239.138 


13 - 


-12 


80 


8.1(-4) 


0.73 (0.26)* 


1.72(0.52)* 


4.60(1.15) 


330.843 


18 6 - 


-17 6 


408 


1.9(-3) 


-0.31 (0.26) 


4.63 (1.28) 


7.44 (2.20) 


330.913 


18 5 - 


-17 5 


329 


2.0(-3) 


- 


3.59(1.13) 


5.91 (1.87) 


330.970 


I84- 


-17 4 


265 


2.K-3) 


- 


2.81 (0.78) 


6.70 (2.07) 


331.014 


I83- 


-17 3 


215 


2.1(-3) 


- 


4.66 (1.30) 


8.26 (2.38) 


331.046 


18 2 - 


-17, 


180 


2.2(-3) 


- 


4.51(1.31) 


10.24 (2.87) 


331.065 


18,- 


-17! 


158 


2.2(-3) 


- 


4.01 (1.09)* 


6.25 (1.96) 


331.072 


18o- 


-17o 


151 


2.2(-3) 


- 


4.01 (1.09)* 


7.86 (2.30) 



238.855 


13 6 - 


-12 6 


337 


9.2(-4) 


- 


<0.22 


<10.23 (2.58): 


238.905 


13 s - 


-12, 


259 


1.0(-3) 


<0.21 


<0.22 


<1.00(0.48): 


238.946 


13 4 - 


-12 4 


195 


l-U-3) 


<0.21 


<0.22 


2.21 (0.70) 


238.978 


13 3 - 


-12 3 


145 


1.1C-3) 


<0.21 


-0.35(0.31) 


1.51 (0.58) 


239.001 


13 2 - 


-12, 


109 


1.2(-3) 


<0.21 


<0.22 


-0.62 (0.35) 


239.015 


13,- 


-12, 


87 


1.2(-3) 


<0.21 


<0.22 


4.61 (0.97)* 


239.020 


13 n - 


-12„ 


80 


1.2(-3) 


<0.21 


<0.22 


4.61 (0.97)* 


330.806 


18 4 - 


-17 4 


265 


3.0(-3) 


<0.31 


<0.32 


<8.08 (2.52): 


330.851 


I83- 


-17 3 


215 


3.1(-3) 


<0.31 


<2.66 (0.53): 


<4.12(1.23): 


330.882 


18 2 - 


-17, 


179 


3.1(-3) 


<0.31 


<0.32 


-0.69 (0.49) 


330.901 


18,- 


-17, 


158 


3.2(-3) 


<0.31 


< 1.49 (0.30): 


2.37 (0.93) 


330.908 


18o- 


-17o 


151 


3.2(-3) 


<0.31 


< 1.49 (0.30): 


3.48(1.12) 



The notation a{-b) stands for a x 
< means upper limit 

- means uncertain detection, S/N 

- means frequency not observed 
: means blended line 

* calculated flux (see Eq. 4) 



<2 



Table A.7: Observed line fluxes J T MB dV (K km s" 1 ) for C 2 H 5 CN. 



Frequency 
[GHz] 



Transition 



F 

'-'Up 

[K] 



A 
[s- 1 ] 



Sources 

IRAS20 126+4 104 IRAS1 



1732 G31.41+0.31 



C,H 5 CN 



218.390 


24 3 , 21 - 


- 23 3-2 o 


140 


8.7(-4) 


<0.18 


- 


-1.03(0.45) 


219.903 


12 3 .io - 


-11 2 ,9 


44 


3.0(-5) 


<0.28 


<0.24 


-0.18(0.13) 


225.236 


25 4 . 2 , - 


- 24 4 , 20 


158 


9.4(-4) 


<0.00 


-0.53 (0.36) 


2.53 (0.85) 


225.307 


12 3 .9 - 


ll2,io 


44 


3.2(-5) 


<0.19 


<0.25 


<1.46(0.57): 


225.317 


23 2-22 - 


- 22 12 i 


122 


4.8(-5) 


<0.19 


<0.25 


<1.79(0.65): 


240.699 


I69.7/8 


- 178,10/9 


148 


7.9(-6) 


<0.20 


<0.48(0.10): 


<0.89 (0.46): 


240.861 


28,. 2 g - 


" 27 ,27 


169 


1.0(-4) 


<0.20 


<0.27 (0.05): 


-0.78 (0.37) 


338.143 


37 3 , 34 - 


- 36 3 , 33 


317 


3.3(-3) 


- 


- 


5.74 (2.44) 


339.895 


39 2 , 38 - 


- 38, 3 7 


334 


3.3(-3) 


<0.15 


-0.38 (0.68) 


3.45(1.05) 


339.968 


38 2 . 3 6 - 


- 37 2 , 35 


327 


3.3(-3) 


<0.15 


-0.50 (0.32) 


2.20 (0.69) 


340.149 


391.38 - 


- 38i_37 


334 


3.3(-3) 


<0.15 


-0.26(0.17) 


2.20 (0.82) 


352.992 


212,19 " 


- 20,,2o 


105 


1.4(-5) 


<0.13 


<0.24 


<0.19 


353.089 


23 3 , 2 o - 


- 22 2 ,2i 


129 


5.8(-5) 


<0.16 


<0.24 


-0.58 (0.29) 


354.477 


40 3 , 38 - 


- 39 3 , 37 


361 


3.8(-3) 


<0.18 


<3.38 (0.68): 


1.33 (0.65) 



The notation a(-b) stands for a x 10 _i 
< means upper limit 

- means uncertain detection, S/N < 2 

- means frequency not observed 
: means blended line 
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Table A.8: Observed line fluxes j T MB dV (K km s _1 ) for HCOOCH3. 



Frequency 


Transition 


p 

-* — Lip 


A 


Sources 






[GHz] 




[K] 


[s- 1 ] 


IRAS20126+4104 


IRAS18089-1732 G31.41+0.31 


HCOOCH, 


218.281 


173,14,2 - 163,13,2 


100 


1.5 (-4) 


<0.18 


- 


1.38(0.42) 


218.298 


173,i4,o - 163,13,0 


100 


1.5 (-4) 


<0.18 


- 


1.00 (0.430) 


219.584 


l°13,5/6,3 - 17l3,4/5,3 


401 


7.7(-5) 


<1.33 


- 


-0.17(0.15) 


219.592 


289,19,2 - 288,20,2 


295 


1.6(-5) 


<1.33 


- 


-0.30 (0.22) 


219.623 


18l2,6/7,3 - 1 '12,5/6,3 


384 


8.9(-5) 


<1.33 


- 


-0.66 (0.32) 


219.642 


18l3,6,4 - 17l3,5,4 


401 


7.7(-5) 


<1.33 


<0.21 


-0.55 (0.43) 


219.696 


18n,8/7,3 - 17n,7/6,3 


369 


1.0(-4) 


<0.28 


<0.21 


0.69 (0.25) 


219.705 


184,15,3 - 174,14,3 


300 


1.5 (-4) 


<0.28 


<0.21 


0.77 (0.27) 


219.764 


189,9,5 - 179,8,5 


342 


1.2(-4) 


<0.28 


<0.21 


0.61 (0.24) 


219.822 


18l0,9/8,3 - 17|0.8/7.3 


355 


1.1 (-4) 


<0.28 


<0.21 


<1. 10 (0.55): 


222.149 


186,12,3 - 176,11,3 


312 


1.5 (-4) 


<0.22 


1.25(0.25): 


<3.05 (0.83): 


222.177 


186,12,5 - 176,11,5 


312 


1.5 (-4) 


<0.18 


<0.20 


-1.14(0.67) 


222.421 


188,10,2 - 178,9,2 


144 


1.3(-4) 


<0.22 


< 1.32 (0.26) 


1.62(0.45) 


222.438 


188,11,0 - 178,10,0 


144 


1.304) 


<0.22 


0.35 (0.16)* 


1.07 (0.39)* 


222.440 


188,10,0 - 17s.9,0 


144 


1.304) 


<0.22 


0.35 (0.16)* 


1.07 (0.39)* 


222.442 


188,11,1 - 178,10,1 


144 


1.304) 


<0.22 


0.35 (0.16)* 


1.07 (0.39)* 


225.372 


2021,9,3 - 1921,8,3 


307 


1.704) 


<0.19 


<0.31 


1.38(0.56) 


225.449 


20n,c),3 ^ 19n,8,3 


307 


1.704) 


<0.19 


<0.18 


<1.30(0.55): 


238.927 


203,18,1 - 192,17,2 


128 


2.105) 


<0.21 


<0.22 


-0.40 (0.23) 


238.933 


203,18,0 ~ 192,17,0 


128 


2-U-5) 


<0.21 


<0.22 


-0.47 (0.20) 


238.947 


193,16,5 - 183,15,5 


309 


2.0(-4) 


<0.21 


-0.29 (0.26) 


<0.88 (0.38): 


239.111 


76,2,4 - 65,2,4 


227 


2.705) 


<0.21 


<0.22 


<0.15 (0.17) 


241.059 


303,27,0 - 303,28,0 


281 


5.8(-6) 


- 


<0.23 


-0.33 (0.31) 


241.068 


303,27,0 - 302,28,0 


281 


1.105) 


- 


<0.23 


-0.76 (0.53) 


330.941 


302/1,29,4/5 - 292/1,28,4/5 


443 


5.504) 


<0.31 


-0.73 (0.38) 


-0.75 (0.83) 


331.021 


255,21,0 - 244,20,0 


210 


3.505) 


<0.34 


<0.40 


<8.29(1.66): 


331.036 


255.21,1 - 24 4 ,20,2 


210 


3.505) 


<0.32 


<0.34 


< 10.01 (3.05): 


331.121 


27i2,16,4 ~26i2,15,4 


505 


4.5(4) 


<0.29 


- 


<0.74 (0.69) 


331.149 


284,25,1 - 27 4 ,24.l 


248 


5.304) 


<0.29 


0.69(0.31) 


2.60(1.39) 


331.160 


284,25,0 - 274,24,0 


248 


5.304) 


<0.29 


<0.99 (0.42) 


<2.24(1.15): 


331.161 


27n, n/16,3 - 26[i, 16/15.3 


490 


4.6(4) 


<0.29 


<0.99 (0.42) 


<2.24 (1.15): 


338.338 


27s,19,2 - 268,18,2 


267 


54(4) 


<0.26 


<0.39 


6.26 (2.37) 


338.356 


27s, 19,0 - 268,18,0 


267 


54(4) 


<0.23 


<0.44 


5.09 (2.08) 


338.393 


285.24,3 - 275,23,3 


443 


5.7(4) 


<0.23 


<0.43 


<3.80(2.05): 


338.396 


277,21,1 - 267.20,1 


258 


5.5(4) 


<0.23 


<0.45 


<3.80 (2.05): 


338.414 


277,21,0 - 267,20,0 


258 


5.5(4) 


<8.06 (2.04): 


<0.45 


<43.32 (9.59): 


339.882 


293,26,3 - 283,25,3 


450 


5.8(4) 


<0.22 


<0.79(0.16) 


-0.73 (0.33) 


340.044 


294.26,4 - 284.25,4 


450 


5.8(4) 


<0.15 


<0.12 


<3.75 (0.75): 


340.115 


277,20,5 - 26 7 19,5 


444 


5.5(4) 


<0.15 


-0.30(0.31) 


0.77 (0.32) 


345.068 


28l4,l4,2 — 27(4,13,2 


370 


4.7(4) 


<0.26 


0.28(0.11) 


<3.93 (1.62): 


345.069 


2814,15/14,0 - 27i4,i4/i3,o 


370 


4.7(4) 


<0.26 


0.56 (0.22) 


<3.93(1.62): 


345.073 


166,11,0 - 155,10,0 


104 


3.9(-5) 


<0.26 


<0.19 


<5.62(1.29): 


345.085 


19l3,6,2 - 19l2,7,2 


224 


3.5(-5) 


- 


-0.37 (0.20) 


<1.73 (0.35): 


345.091 


2814,15,1 -27 14 ,i 4 j 


370 


4.7(4) 


<0.26 


<0.50 (0.32) 


<4. 14 (1.70): 


345.148 


286,23,3 - 27 6 ,22,3 


452 


6.0(4) 


<0.26 


-0.55 (0.34) 


-0.59 (0.54) 


345.163 


118,3,5 - 107,3,5 


269 


7.1(-5) 


<0.26 


<0.19 


-0.50 (0.47) 


345.230 


18l3,5,2 - 18l2,6,2 


213 


3.2(-5) 


<0.25 


-0.59 (0.40) 


<8.18(3.04): 


345.242 


18l3,6,l - 18l2,7,l 


213 


3.2(-5) 


<0.25 


<0.85 (0.83) 


<8.18(3.04): 


345.248 


28l0,19,4 -27io,i8,4 


493 


5.5(4) 


<0.26 


<0.19 


<0.34 


352.817 


33o/i/o/i,33,3 - 32i/[/o/o,32,3 


479 


1.1(4) 


<0.18 


-0.38 (0.28) 


<2.68 (0.76): 


352.841 


33l/o,33,4/5 - 32i/o,32,4/5 


479 


7.7(4) 


<0.18 


<0.50 (0.42) 


1.72(0.47) 


352.912 


312,29,2 - 303.28,1 


286 


8.5(-5) 


<0.18 


<0.26 


0.30 (0.08) 


352.918 


313,29,1 - 303,28,1 


286 


6.6(4) 


<0.18 


<1. 10 (0.36) 


2.32 (0.62)* 


352.922 


312,29,2 - 302,28,2 


286 


6.6(4) 


<0.18 


<1. 10 (0.36) 


2.32 (0.62)* 


352.926 


313,29,0 - 303.28,0 


286 


6.6(4) 


<0.18 


<1. 10 (0.36) 


2.32 (0.62)* 


352.930 


312.29,0 ~ 302.28,0 


286 


6.6(4) 


<0.18 


<1. 10 (0.36) 


2.32 (0.62)* 



The notation a(-b) stands for a x 10~* 
< means upper limit 

- means uncertain detection, S/N < 2 

- means frequency not observed 
: means blended line 

* calculated flux (see Eq. 4) 
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Table A.9: Observed line fluxes j T MB dV (K km s" 1 ) for HCOOCH3 continued. 



Frequency 


Transition 


^ — Lip 


A 


Sources 






[GHz] 




LKJ 


[s- 1 ] 


IRAS20 126+4 104 


IRAS18089-1732 


G31.41+0.31 


HCOOCH3 


354.427 


29i6,13/14,3 - 28i6,12/13,3 


614 


4.7(-4) 


<0.18 


<0.23 


-1.04(0.64) 


354.477 


29l5_l5,4 -28i5j4_4 


593 


5.0(4) 


<0.13 


<0.23 


<0.40 (0.58): 


354.574 


128.5/4,3 - 117,4/5,3 


276 


7.405) 


<0.18 


<0.23 


< 1.44 (0.78): 


354.605 


29l5_l5/l4,3 - 281544/13,3 


593 


5.0(4) 


<0.18 


<0.23 


1.96(0.46) 


354.608 


33l/o,33,l/2 - 32]/o,32,l/2 


293 


7.2(4) 


<0.18 


1.11(0.42)* 


3.14(0.49)* 


354.608 


33o/l,33,2/l - 32i/o,32,l/2 


293 


6.8(4) 


<0.18 


1.05 (0.40)* 


2.75 (0.47)* 


354.629 


287.21,5 - 27 7 . 2 0,5 


461 


6.4(4) 


<0.18 


<0.23 


-1.20(0.68) 


354.742 


128,5,1 - H7,5.1 


88 


7.3(-5) 


<0.18 


<0.23 


1.19(0.44) 


354.759 


128,4,2 - H7,4,2 


88 


7.3(-5) 


<0.16 


<0.17 


<3.23 (0.65): 


354.806 


128.5/4,0 - 117,4/5,0 


88 


7.3(-5) 


<0.18 


<0.23 


3.90(1.18) 


354.839 


29i4_i5/i6_3 - 2814,14/153 


574 


5.2(4) 


<0.15 


<0.36 (0.46): 


<4.91 (1.40): 


364.297 


33l/2,32,2/l - 321,31,2 


308 


7.3(4) 


<0.24 


<6.81 (1.36): 


<3.67 (1.16): 


364.302 


332/1,32,0 - 32 2 /l,3l,o 


308 


7.3(4) 


0.80(0.91) 


<6.81 (1.36): 


<5.00(1.42): 



The notation a(-b) stands for a x 10~* 
< means upper limit 

- means uncertain detection, S/N < 2 

- means frequency not observed 
: means blended line 

* calculated flux (see Eq. 4) 



Table A.10: Observed line fluxes fT MB dV (K km s _1 ) for CH 3 OCH 3 . 



Frequency 
[GHz] 



Transition 



[K] 



A 
[s- 1 ] 



Sources 
IRAS20126+4104 



IRAS 18089-1732 G31.41+0.31 



CH3OCH3 



The notation a(-b) stands for a x 10~* 
< means upper limit 

- means uncertain detection, S/N < 2 

- means frequency not observed 
: means blended line 

* calculated flux (see Eq. 4) 



218.490 


233,21,2/3 - 232,22,2/3 


264 


3.4(-5) 


<0.18 


- 


<1. 19 (0.75) 




218.492 


233,21,1 - 232,22,1 


264 


3.4(-5) 


<0.18 


- 


<1. 19 (0.75) 




218.495 


233,21,0 - 232,22,0 


264 


3.4(-5) 


<0.18 


- 


<1. 19 (0.75) 




222.239 


4-3,2,2 - 32,1,2 


22 


2.8(-5) 


<0.22 


0.82 (0.56) 


-0.93 (0.43) 


222.248 


43,2.3/1 - 32,1,3/1 


22 


4.2(-5) 


<0.22 


-0.43 (0.20) 


<1. 15 (0.35): 


222.255 


43,2,0 - 32,1.0 


22 


4.9(-5) 


<0.22 


-0.62 (0.21) 


1.11 (0.98) 


222.259 


43,1,1 - 32,1,1 


22 


1.6(-5) 


<0.22 


-0.30 (0.21)* 


2.37 (0.74)* 


222.260 


4-3,1,2 - 32,1,2 


22 


2.1(-5) 


<0.22 


-0.11(0.21)* 


0.83 (0.26)* 


222.326 


253.23,0 - 244.20,0 


308 


7.9(-6) 


<0.22 


<0.27 


-0.42(0.19)* 


222.327 


253.22,1 - 24 4 .20,1 


308 


7.9(-6) 


<0.22 


<0.27 


-1.12(0.52)* 


222.329 


253,22/23.2/3 - 244,20/20,2/3 


308 


7.9(-6) 


<0.22 


<0.27 


-0.42 (0.21)* 


222.414 


43,2,2 - 3i,2,2 


22 


2.1(-5) 


<0.22 


-0.20(0.10) 


<0.73(0.15) 




222.423 


43,2.1 - 32,2.1 


22 


1.6(-5) 


<0.22 


-0.56 (0.30) 


<2.09 (0.42) 




222.427 


43,1.3 - 32,2.3 


22 


4.9(-5) 


<0.22 


-0.23(0.12) 


<0.85(0.11) 




222.434 


43,1,1/0 - 32,2,1/0 


22 


4.2(-5) 


<0.22 


1.68(0.63) 


<1.86 (0.52) 




222.435 


43,1,2 - 32,2,2 


22 


2.8(-5) 


<0.22 


-0.22 (0.08) 


<0.34(0.10) 




225.202 


244,21,2/3 - 243,22,2/3 


296 


4.7(-5) 


<0.19 


<0.30 


0.55 (0.18)* 


225.204 


244.21,1 - 24 3 ,22,1 


296 


4.7(-5) 


<0.19 


-0.15(0.14)* 


0.87 (0.28)* 


225.205 


244,21,0 - 243,22,0 


296 


4.7(-5) 


<0.19 


-0.09 (0.08)* 


0.55 (0.18)* 


238.975 


295,25,2/3/1/0 - 286,22,2/3/1/0 


432 


1.6(-5) 


<0.21 


<0.22 


< 1.64 (0.64) 




239.020 


245.19,5/3 - 244,20,5/3 


309 


6.1(-5) 


<0.21 


<0.63(0.12): 


<0.68(0.14) 




239.021 


245.19,1 - 24 42 0,1 


309 


6.1(-5) 


<0.21 


< 1.67 (0.33): 


< 1.82 (0.36) 




239.021 


245,19,0 - 244,20,0 


309 


6.1(-5) 


<0.21 


<0.63(0.13): 


<0.68(0.14) 




240.978 


53,3,2 - 42,2.2 


26 


4.5(-5) 


<0.20 


-0.11(0.07)* 


2.16 (0.66)* 


240.983 


53,3,3 - 42,2.3 


26 


5.4(-5) 


<0.20 


-0.06 (0.04)* 


0.19 (0.06)* 


240.985 


53,3,1 _ 4 2 ,2.1 


26 


5.1(-5) 


<0.20 


-0.49 (0.30)* 


1.49 (0.46)* 


240.990 


53,3,0 - 42,2,0 


26 


5.4(-5) 


<0.20 


-0.19(0.12) 


<0.64 (0.28): 
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Table A. 1 1 : Observed line fluxes j T MB dV (K km s _1 ) for CH 2 CO. 



Frequency 
[GHz] 



Transition 



'-'Up 

[K] 



A 
[s" 1 ] 



Sources 
IRAS20126+4104 



IRAS18089-1732 G31.41+0.31 



CH 2 CO 


















222.120 


114,8/7 


- IO4/7/6 


273 


1.1 (-4) 


<0.22 


<0.27 


<0.54(0.10) 




222.198 


llo.11 - 


- lOo.io 


65 


1.2(-4) 


<0.22 


<0.81 (0.26): 


<2.34 (0.75) 




222.201 


113,9/8 


- 103,8/7 


181 


1.1 (-4) 


<0.22 


<0.81 (0.26): 


<2.34 (0.75) 




222.229 


112,10 - 


- 10 2 ,9 


116 


1.2(-4) 


<0.22 


<0.28 (0.06): 


<2.22 (0.32) 




222.315 


112,,- 


10 2 ,8 


116 


1.2(-4) 


<0.22 


-0.24 (0.22) 


-0.91 (0.57) 


363.937! 


182,16 " 


- 17 2 ,15 


218 


5.5(-4) 


<4.65 (1.08): 


<0.22 


<20.29(3.31): 



The notation a(-b) stands for a x 10~' 

< means upper limit 

~ means uncertain detection, S/N <2 

- means frequency not observed 

: means blended line 

! frequency from CDMS 



Table A.12: Observed line fluxes / T MB dV (K km s -1 ) for CH3CHO. 



Frequency 
[GHz] 



CH3CHO 



Transition 



[K] 



A 
[s" 1 ] 



Sources 
IRAS20126+4104 



IRAS 18089- 1732 G3 1.41 +0.31 



330.822 


53,3,1 - 42,3,1 


34 


1.2(-4) 


<0.31 


-0.71 (0.80) 


-0.97 (0.52) 


331.039 


173,14,3 - 163,13,3 


367 


1.2(-3) 


<1. 18 (0.88): 


<0.42 (0.29): 


<2.50(1.47): 


354.458 


182,16,3 - 172,15,3 


375 


1.5(-3) 


<0.18 


<3.11 (1.11): 


<3.25(1.26): 


354.525 


19o,19,5 - 1 80, 18,5 


377 


1.5(-3) 


<0.17 


<0.19 


<0.34 


354.813 


182,16,2 - 172,15,2 


170 


1.5(-3) 


<0.18 


<0.23 


-1.54(0.92) 


354.844 


182,16,0 - 172,15,0 


170 


1.5(-3) 


<0.18: 


<1. 13 (0.55): 


<4.92(1.34): 



The notation a(-b) stands for a x 10 
< means upper limit 

- means uncertain detection, S/N < 2 

- means frequency not observed 
: means blended line 



Table A. 13: Observed line fluxes J TMBdV (K km s ') for HCOOH and its isotopic species. 



Frequency 


Transition 


p 


A 


Sources 






[GHz] 




[K] 


[s- 1 ] 


IRAS20126+4104 


IRAS18089-1732 


G3 1.4 1+0.31 


HCOOH 


222.110 


72,6 _ 7l,7 


43 


1.9(-6) 


<0.22 


<0.27 


<0.24 


225.086 


10 4 ,7-9 4 ,6 


110 


1.0(-4) 


<0.18 


<0.18 


<1.97 (1.70): 


225.091 


10 4 ,6-9 4 ,5 


110 


1.0(-4) 


<0.18 


<0.18 


<1.97 (1.70): 


225.238 


10 3 ,8-9 3 ,7 


88 


1.1 (-4) 


<0.19 


<0.39 (0.24): 


<2.59 (0.78): 


330.931 


4 3 ,i - 4 2 ,2 


39 


5.1 (-6) 


<0.31 


<0.32 


-0.47 (0.63) 


331.145 


33,0 - 32,1 


35 


3.6(-6) 


<0.29 


<0.29: 


<2.77 (1.80): 


338.109 


15l,15 - 14o,14 


127 


l.l(-5) 


- 


- 


<5.11 (3.10): 


338.202 


153,13 - 143,12 


158 


4.1 (-4) 


- 


- 


<13.55 (6.28): 


345.031 


16o,16 - 15o,15 


143 


4.5(-4) 


<0.26 


-0.30 (0.30) 


-1.18(0.70) 


345.253 


14 3 ,12 - 14 2 ,13 


142 


8.3(-6) 


<0.29 


<0.19: 


<1.59 (1.01): 


354.448 


17o,17 - 16[,16 


161 


1.3(-5) 


<0.18 


<1.70(0.57): 


<6.42(1.60): 



The notation a(-b) stands for a x 10"' 
< means upper limit 

- means uncertain detection, S/N <2 

- means frequency not observed 
: means blended line 
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Table A.14: Observed line fluxes / T MB dV (K km s 1 ) for CH 3 CCH. 



Frequency 


Transition 




A 


Sources 






[GHz] 






[K] 


[s- 1 ] 


IRAS20126+4104 


IRAS18089-1732 


G31.41+0.31 


CH 3 CCH 


222.099 


13 4 - 


-12 4 


190 


3.4(-5) 


<0.22 


<0.65 (0.35): 


<2.79 (0.56): 


222.129 


13a- 


-12 3 


140 


3.6(-5) 


-0.29(0.17) 


1.14(0.43) 


3.15 (0.94) 


222.150 


132- 


-122 


104 


3.7(-5) 


-0.41 (0.19) 


1.25 (0.46) 


<4.25 (0.85): 


222.163 


13i- 


-12i 


82 


3.8(-5) 


0.74 (0.29) 


2.20 (0.68) 


3.14(1.01) 


222.167 


13 - 


-12 


75 


3.8(-5) 


0.93 (0.35) 


2.13(0.70) 


3.27 (0.97) 


239.088 


IV 


-13 6 


346 


3.9(-5) 


<0.21 


<0.22 (0.05): 


<2.40 (0.48): 


239.138 


IV 


-13 5 


267 


4.1(-5) 


<0.21 


<2.38 (0.48): 


<6.39(1.28): 


239.179 


IV 


-13 4 


202 


4.3(-5) 


0.71 (0.14) 


<0.68(0.14): 


-0.39 (0.21) 


239.211 


IV 


-13 3 


151 


4.5(-5) 


-0.40 (0.24) 


-0.75 (0.41) 


2.47 (0.64) 


239.234 


IV 


-132 


115 


4.6(-5) 


0.36(0.16) 


-0.89 (0.89) 


2.06 (0.63) 


239.248 


IV 


-13! 


93 


4.7(-5) 


0.79 (0.23) 


1.62(0.52) 


2.91 (0.74) 


239.252 


IV 


-13o 


86 


4.7(-5) 


0.87 (0.24) 


1.68(0.43) 


3.46(1.07) 



The notation a(-b) stands for a x 10 
< means upper limit 

- means uncertain detection, S/N <2 

- means frequency not observed 
: means blended line 
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Appendix B: Rotation diagrams 

RTD diagrams for H 2 CO, CH 3 OH, C 2 H 5 OH, HNCO, NH 2 CHO, 
CH 3 CN, C2H5CN, HCOOCH3, CH3OCH3, CH 2 CO and 
CH3CCH. Lines from different frequency bands are corrected 
for differential beam dilution assuming a Rt=M)k source size for 
the warm species and 14" source size for the cold species, as in- 
dicated in the plots. Optically thin, unblended lines with S/N>2 
are marked with filled circles. Lines with high S/N ratio (<2) 
are marked with diamonds, those included in the fit with filled 
diamonds. Optically thick lines and blended lines are marked 
with open squares and triangles, respectively. Upper limits are 
marked with arrows. Upper limits used to constrain the fit are 
marked with filled stars. The error bars are calculated using Eq. 
A.2. 
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24 



22 



T rot = 157 +/- 44 K 



(1.28 +/- 0.45) e+15 cm-2 
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E up (K) 



28 




IRAS1 8089- 1732, H 2 C0, 2.4". 


26 




^~r^^^ : 


24 




^^-^ 




■ T rot = 


204 +/- 82 K 


?? 


^beom 


= (5.14 +/- 2.00) e+14 cm-2 . 



100 200 300 400 500 
E up (K) 



28 



IRAS20 126+4 104, H 2 C0, 2.2' 




100 200 300 400 500 
E up (K) 



Fig. B. 1 : RTD fits for H 2 CO. All lines included into the fit belong 
to para-H 2 CO. 



Isokoski et al.: m09bnl0. 




(2.32 +/- 0.34) e+16 cm-2 
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E UP (K) 





IRAS201 26+41 04, CH 3 0H, 2.2" '. 


-t 


f : 






I. 

v 4- 4-. 

4 i * * i. 


- T rot = 


119 +/- 17 K^ 


\ 4. . 


'^beom 


= ( 7.09 +/- 2.76) e+>KQrn-i 



200 400 600 800 

E up (K) 

Fig. B.2: RTD fits for CH 3 OH. Lines with £ up <100 K are con- 
sidered contaminated with the cold CH3OH. Lines with S/N>1 
are included in the fit. 
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120 +/- 15 K : 
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= (2.87 +/- 0.67) e+15 cm-2 : 



100 200 300 400 500 
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30 
28 



IRAS1 8089- 1732, C 2 H 5 0H, 2.4" 




= (1.06 +/- 0.66) e+15 cm-2 



100 200 300 400 500 
E up (K) 

Fig.B.3: RTD fits for C 2 H 5 OH. No C 2 H 5 OH was detected in 
IRAS20126+4104. ForIRAS18089 lines with S/N>1 have been 
used. 
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94 +/- 50 K J 


^beom 
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Fig.B.4: RTD fits for HNCO. No HNCO was detected in Fig. B.5: RTD fits for NH 2 CHO. No NH 2 CHO was detected in 



IRAS20126+4104. 



IRAS20126+4104. 



G31.41+0.31, CH 3 CN, 2.0" 
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Fig.B.7: RTD fits for C 2 H 5 CN. No C 2 H 5 CN was detected in 
IRAS20126+4104. For IRAS18089 lines with S /N < 2 have 
been included in the fit. 
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Fig. B.6: RTD fits for CH 3 CN. 
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Fig. B.8: RTD fits for HCOOCH 3 . No HCOOCH 3 was detected Fi B g. RTD fits for C H 3 OCH 3 . No CH 3 OCH 3 was detected in 
in IRAS20126+4104. IRAS20126+4104. 
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Fig. B. 10: RTD fits for CH 2 CO. No CH 2 CO was detected in 
IRAS20126+4104. 
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Fig. B. 12: RTD fits for HCOOH. 
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Fig. B.l 1: RTD fits for CH3CHO. 
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Fig. B.13: RTD fits for CH 3 CCH. 
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Appendix C: Weeds model parameters Table C1: Weeds model parameters for IRAS20126+4104 



Species 


N s 


T 

a ex 


o s 


Av 




[cm- 2 ] 


[K] 


["] 


[km s" 1 ] 


CH3OH 


1.5E+17 


300 


2.2 


6 


CH3OH 


4.8E+16 


20 


3.0 


5 


13 CH 3 OH 


- 


- 


- 


- 


H 2 CO 


6.0E+15 


150 


2.2 


6 


H 2 13 CO 


<2.0E+15 


150 


2.2 


6 


CH 3 CN ' 


1.3E+15 


200 


2.2 


6 


CH 3 13 CN 


<2.0E+13 


200 


2.2 


6 


HNCO 


1.0E+15 


200 


2.2 


6 


HN 13 CO 


<1.0E+14 


200 


2.2 


6 


C2H5OH 


<3.0E+15 


100 


2.2 


6 


C 2 H 5 CN 


<5.0E+14 


80 


2.2 


5 


CH3OCH3 


<1.0E+16 


100 


2.2 


6 


HCOOCH3 


<2.0E+15 


200 


2.2 


6 


NH2CHO 


3.0E+14 


300 


2.2 


6 


CH 2 CO 


<5.0E+13 


50 


14.0 


6 


CH3CHO 


<1.0E+14 


50 


8.0 


6 


HCOOH 


<5.0E+13 


40 


8.0 


6 


CH3CCH 


7.0E+14 


35 


14.0 


2 


Additional species 


CH3COCH3 


<0.1E+17 


300 


2.2 


5 


HCCCN 


1E+13 


100 


14 


7 


SO 


0.3E+15 


50 


14 


6 


34 SO 


0.2E+14 


50 


14.0 


6 


so 2 


0.3E+15 


50 


14.0 


6 


33 S02 


0.05E+15 


50 


14.0 


5 


HCN 


0.4E+14 


50 


14.0 


20 


HCN 


4E+13 


50 


14.0 


5 


H 13 CN 


0.1E+13 


50 


14.0 


6 


CN, v = 0, 1 


0.2E+15 


50 


14 


2 



13 CH 3 OH is not available in the JPL database used for Weeds modeling 
1 based on 21" beam spectra. 



Isokoski et al.: m09bnl0, Online Material p 20 
Table C.2: Weeds model parameters for IRAS18089-1732 Table C.3: Weeds model parameters for G31.41+0.31 



Species 


N s 


T 

A ex 


0s 


Av 




[cm" 2 ] 


[K] 


["] 


[kins" 1 ] 


CH 3 OH 


3.5E+17 


300 


2.4 


5 


CH 3 OH 


2.6E+17 


20 


3.0 


5 


13 CH 3 OH 


- 


- 


- 


- 


H 2 CO 


1.5E+16 


150 


2.4 


6 


H 2 13 CO 


3.5E+15 


150 


2.4 


6 


CH 3 CN 


3.5E+15 


200 


2.4 


6 


CH 3 13 CN 


4.0E+14 


200 


2.4 


6 


HNCO 


4.0E+15 


200 


2.4 


6 


HN 13 CO 


7.0E+14 


200 


2.4 


6 


C 2 H 5 OH 


2.5E+16 


150 


2.4 


6 


C 2 H 5 CN 


4.0E+15 


80 


2.4 


5 


CH 3 OCH 3 


1.0E+17 


100 


2.4 


6 


HCOOCH 3 


3.0E+16 


200 


2.4 


6 


NH 2 CHO 


5.0E+14 


100 


2.4 


6 


CH 2 CO 


1.5E+14 


50 


14.0 


5 


CH 3 CHO 


<1.0E+14 


50 


8.0 


6 


HCOOH 


<5.0E+13 


40 


8.0 


6 


CH 3 CCH 


1.8E+15 


40 


14.0 


3.5 



Additional species 



CH 3 COCH 3 


0.5E+17 


300 


2.4 


5 


HCCCN 


3E+13 


150 


14 


6 


SO 


1E+15 


50 


14 


6 


34 SO 


0.5E+14 


50 


14.0 


6 


so 2 


0.5E+15 


50 


14.0 


6 


33 so 2 


0.5E+15 


50 


14.0 


5 


HCN 


2E+14 


50 


14.0 


7 


HCN 


-5E+13 


50 


14.0 


4 


H 13 CN 


1.5E+13 


50 


14.0 


6 


CN, v = 0, 1 


0.4E+15 


50 


14 


5 


NH 2 CN 


0.1E+14 


50 


14 


7 


OC 34 S 


0.5E+15 


50 


14.0 


7 


CP 


0.5E+15 


50 


14.0 


4 


c 34 s 


1.4E+14 


50 


14.0 


5 


HCOCH 2 OH 


0.05E+17 


300 


2.0 


6 



13 CH 3 OH is not available in the JPL database used for Weeds modeling. 



Species 


N s 


T 

x ex 


0s 


Av 




[cm" 2 ] 


[K] 


["] 


[km s~'] 


CH 3 OH 


2.0E+18 


300 


2.0 


7 


CH 3 OH 


2.7E+17 


20 


5.0 


7 


13 CH 3 OH 


- 


- 


- 


- 


H 2 CO 


6.0E+16 


150 


2.0 


6 


H 2 13 CO 


7.0E+15 


150 


2.0 


6 


CH 3 CN ' 


2.0E+16 


300 


2.0 


7 


CH 3 13 CN 


2.0E+15 


300 


2.0 


7 


HNCO 


6.0E+16 


200 


2.0 


7 


HN 13 CO 


2.0E+15 


200 


2.0 


7 


C 2 H 5 OH 


2.0E+17 


100 


2.0 


7 


C 2 H 5 CN 


2.0E+16 


80 


2.0 


5 


CH 3 OCH 3 


1.0E+18 


100 


2.0 


5 


HCOOCH 3 2 


1.0E+18 


300 


2.0 


7 


NHjCHO 


4.0E+15 


300 


2.0 


7 


CH 2 CO 


6.5E+14 


50 


14.0 


6 


CH 3 CHO 


1.0E+15 


50 


8.0 


6 


HCOOH 


0.2E+16 


40 


8.0 


6 


CH 3 CCH 


2.8E+15 


60 


14.0 


4 



Additional species 








CH 3 COCH 3 


1E+17 


100 


2.0 


7 


HCCCN 


7E+13 


150 


14 


7 


SO 


7E+14 


50 


14 


8 


34 SO 


1E+14 


50 


14.0 


6 


S0 2 


0.1E+16 


50 


14.0 


7 


33 so 2 


1E+15 


50 


14.0 


6 


HCN 


2E+14 


50 


14.0 


12 


HCN 


-6E+13 


50 


14.0 


6 


H 13 CN 


2E+13 


50 


14.0 


6 


CN, v = 0, 1 


0.8E+15 


50 


14 


6 


NH 2 CN 


0.3E+14 


50 


14 


7 


o 13 cs 


5E+14 


50 


14 


5 


OC 34 S 


1E+15 


50 


14.0 


7 


CP 


1E+15 


50 


14.0 


6 


c 34 s 


2E+14 


50 


14.0 


7 


HCOCH,OH 


0.2E+17 


300 


2.0 


7 


s 3 


1.5E+16 


15 


14 


6 


CH 3 NH 2 


1E+17 


150 


2.0 


7 



13 CH 3 OH is not available in the JPL database used for Weeds modeling. 

1 based on 21" beam spectra. 

2 based on 21" beam spectra, 1.5E+17 in the 14" beam spectra. 
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Appendix D: Additional detections 

Several lines of other species were found in the observed fre- 
quency ranges, particularly in the line-rich G31.41+0.31. Table 
D.l lists additional detections with column densities obtained 
from the Weeds analysis with the best T ex when this could be de- 
rived. Several transitions of acetone, CH3COCH3, are detected. 
In G31, the best agreement is obtained at 100 K, whereas for 
IRAS 18089, this is r ex =300 K. No CH 3 COCH 3 was detected in 
IRAS20126, and the tabulated value is an upper limit at 300 K. 
Figure D.l shows the strongest acetone lines in G31 together 
with the Weeds model. 

Glycoaldehyde, HCOCH2OH, has previously been detected 
in G31 (Beltran et al. 2009), has several lines in the covered 
ranges. All lines are however blended with other transitions, and 
only upper limit could therefore be derived reliably. The tabu- 
lated column density is for a temperature of 300 K, constrained 
by non-detection of lines with low E up . 

Two transitions of HC3N are detected at 218.325 GHz 
(£ up 131.0 K) (J = 24 -> 23) and 354.697 GHz (£ up 340.5 
K) {J - 39 — > 38). The Weeds model on the two lines, well 
separated in E up , gives a beam averaged column density of 
7 x 10 13 cm -2 and a temperature of 150 K for G31. The HC3N 
emission cannot be matched with emission contained within a 
2.0" volume as the low-Euptransition becomes optically thick. 
For IRAS20126, the tabulated column density is at 100 K, con- 
strained by the two transitions. Only the high E up transition was 
covered for IRAS 18089 giving a column density of 3x 10 13 cm -2 
at 150 K and 6 x 10 13 cm -2 at 100 K. 

Several CH3NH2 transitions are observed with a T ex of 
150 K. Blends with other species make the rotation tempera- 
ture and column density inaccurate, however. Several unidenti- 
fied lines coincide with NH2CH2CH2OH, but only upper limits 
of l.OxlO 16 cm -2 at 100-300 K could be derived reliably using 
Weeds. 

Several weak lines of NH2CN are detected, with a T ex of 
50 K derived from the Weeds model, but line blends make the 
rotation temperature and column density inaccurate. For CN, 
SO, 34 SO, 13 CS, OC 34 S, S0 2 , 33 S0 2 , C 34 S, and CP not enough 
lines were observed to derive r ex values, and tabulated column 
densities are assuming a temperature of 50 K. 

One line of HCN, 7=4-3, was observed. HCN emission 
in G31 is composed of a broad (12 km s _1 ) and a narrow (6 
km s _1 ) component (redshifted by 2.5 km s _1 ) causing self- 
absorption. For IRAS 18089, the broad component has a width 
of 7 km s and narrow 4 km s -1 (redshifted by 1 km s -1 ). For 
IRAS20126 the HCN emission has a very broad 20 km s com- 
ponent blueshifted by 2 km s and a narrow (5 km s _1 ) compo- 
nent redshifted by 1 km s . The tabulated column densities are 
those derived from the H I3 CN column density. 



CH3COCH3 

(23i,22,0 - 22221,1^ 
(232,22,0 ~ 22 2 ,2i,i) 
(23i,22,o - 22j 21,1) 
(232,22,0 ~ 22i ,21,1) 




CH3COCH3 

(287,21,0 ~ 27 8 20,1 

(28g,2l,o - 27 7 20,1 




340.09 340.10 340.11 

Rest Frequency (GHz) 

Fig. D.l: Beam-averaged spectra of CH3COCH3 (23i/2,22,i~ 
22i/2,2i,i), (23i/2,22,o-22i/2,2i,i), (287/8,2 i,o-27g/ 7j 2o,o) and 
(287/8 j 2i,o-27g/7 j 2o,i) towards G31.41+0.31 hot core. Dashed 
lines indicate Gaussian fits to the lines and red solid line show 
the Weeds model on the CH3COCH3 lines at 100 K. 
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Table D.l: Additional detections. Column densities given are source averaged for species with r ex >100 K and beam averaged for 
species with r ex <100K (except for HCCCN, the emission of which is warm and arises from extended volume). The tabulated values 
are those obtained from Weeds analysis with the best T ex . 



Species 



IRAS20126+4104 IRAS 18089- 1732 G31.41+0.31 



Source-averaged column densities 



CH3COCH3 


<0.1E+17 


0.5E+17 


1.0E+17 


HCOCH20H 


- 


<0.5E+16 


<0.2E+17 


CH 3 NH 2 


- 


- 


2.4E+14 


NH2CH2CH2OH 


- 


- 


<1.0E+16 



Beam-averaged column densities 



NH 2 CN 


- 


<0.1E+14 


7.0E+14 


OC-13-S 


- 


- 


5.0E+14 


OC-34-S 


- 


0.5E+15 


1.0E+15 


HCCCN 


1.0E+13 


3.0E+13 


7.0E+13 


S3 


- 


- 


1.5E+16 


CN, v = 0, 1 


0.2E+15 


0.4E+15 


0.8E+15 


S0 2 


0.3E+15 


0.5E+15 


0.1E+16 


SO 


0.3E+15 


1.0E+15 


7.0E+14 


CS-34 


- 


1.4E+14 


2.4E+14 


CP 


- 


- 


1.0E+15 


S-33-0 2 


<0.5E+14 


0.5E+15 


1.0E+15 


S-34-O 


0.2E+14 


0.5E+14 


1.0E+14 


HCN 


8E+13 


2.0E+14 


2.0E+14 


HC-13-N 


0.1E+13 


1.5E+13 


2.0E+13 



